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ABSTRACT 


This  thesis  is  concerned  with  the  time 
evolution  of  a  large  system  consisting  of  several 
kinds  of  particles  that  can  be  transformed  into 
each  other  by  an  external  field.  Several  kinds  of 
systems  are  studied,  and  in  each  case  an  exact 
equation  describing  the  time  evolution  of  the  system 
is  obtained. 

Chapter  One  is  an  introduction  to  the  problem. 

It  is  a  recapitulation  of  some  work  done  by  H.J. 

Kreuzer  and  K.  Nakamura  [1] .  They  considered  two 
kinds  of  Bose-Einstein  particles  a  and  b,  which  can 
be  transformed  into  each  other  via  a  constant  external 
field.  We  extend  this  to  include  Fermi-Dirac  particles. 

Chapter  Two  generalizes  the  work  of  Chapter 
One  for  bosons  to  a  three-channel  model  (i.e.  with 
three  kinds  of  particles  a,  b,  and  c) .  The  resulting 
time  evolution  equation  is  then  studied  and  illustrated 
by  various  examples.  The  extension  to  an  n-channel 
system  of  arbitrary  topology  is  easily  made.  A  time 
evolution  equation  is  derived  but  not  studied. 

The  zero  temperature  limit  for  fermions  is 
obtained  in  Chapter  Three.  Kelvin's  stationary  phase 
argument  is  used  to  study  the  asymptotic  time  evolu¬ 
tion  . 
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In  Chapter  Four  we  study  a  two-channel  system 
initially  in  equilibrium  but  with  the  a  particles  at 
temperature  T  and  the  b  particles  at  temperature  T,  . 

^  Ij 

Finally  in  Chapter  Five  a  high  temperature 
expansion  is  made  for  weakly  degenerate  two-channel 
boson  and  fermion  systems. 
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CHAPTER  I 


THE  TWO  CHANNEL  SYSTEM 

We  consider  two  kinds  of  particles  a  and  b 
(either  both  bosons  or  both  fermions)  which  can  be 
transformed  into  each  other  via  an  external  field. 


The  hamiltonian 

The  hamiltonian  for  this  sytem  in  coordinate 
representation  is 


J4  ~ 


2m, 


V^-V^dx  + 


-> 

Vip 2  •  Vij^dx 


+  e  (t) 


dx  V  (x)  [ip^2  + 


(1.1) 


where  (x,t),  i=  1,2  are  either  both  Bose  or 

both  Fermi  field  operators  and 


6  (t) 


1 


t  <  0 
t  >  0 


Fourier  decomposing  the  field  operators,  we  get 


ip1  (x ,  t )  =  —  l  a£  e 
/V 

k 


•  T*" 

i  k .  x 


(1.2) 


\{j?  (x ,  t )  =  —  X  b^r 


i  k.x 


/V 


1 


' 


■ 


2 


NOTE  Henceforth,  the  wavenumber  k  will  simply  be 
written  k,  the  dot  product  k.x,  etc.  Also 


I  =  l 

k  (kx,ky,kz) 


where 


27rn . 
k  =  _ i 

1  ^73  ' 


l  =  x ,  y ,  z  ;  n^  =  0, 1, . . . 


OO 


Substituting  (1.2)  into  (1.1),  we  get  the  hamiltonian 
in  creation-  and  annihilation-operator  representation 


l  Ek1)akak+|  ek2>bkbk+  6(t)kJk,Vkk'  (akbk-+b^ak,) 


(1.3) 


where 


(i) 


ft2k2 

2m . 


V 


kk'  V 


dx  V  (x)  e 


i (k-k ' ) -x 


Note  that  if  V (x)  =  Vq  (a  constant  interaction)  then 
Vkk'  = 


o  kk 


Diagonalization  and  exact  solution 


In  matrix  notation,  the  hamiltonian  equation 


is 


V 


^  =  (aklbk) 


V 


t 


(2) 


(1.4) 


. 


where  E 


(i) 


and  V  are  submatrices  with  elements 

t,Lt. 


Ekk'  =  £kl)(5kk'  and  Vkk '  resPectively ;  and  (a^|b^) 

4.  4. 

represents  the  infinite  vector  ( .  .  .  a  '  . . .  |  .  b '  .  .  .  ) 

K  K 

with  all  wavenumbers  k  included.  The  hamiltonian 
equation  (1.4)  can  be  diagonalized  with  a  unitary 
matrix 

r 

A  B 

U  = 

C  D 


That  is 


=  I  (a£|b^)U  X] 


+  t  ak 

(ak  |  3k)  A  (g-) 


U 


'e(1)  V 


e(2) 


u 


-1 


u 


^bk' 


(1.5) 


where 


U 


E(l) 

> 

V 

-1 

/ 

o 

1 — 1 

< 

v. _ 

v+ 

e(2> 

J 

u  = 

0  a(2) 

=  A 


(1.6) 


and  A 

A 


(i) 


are  diagonal  submatrices  with  elements 


(i)  _  ,(i) 


kk 


We  introduce  quasi-particles  by  a 


linear  transformation  of  the  annihilation  operators 


ak  \x  (Akk,ak'  +  Bkk,bk' } 

k ' 


(1.7) 


3k  (Ckk,ak'  +  Dkk,bk' } 


with  inverse  transformations 


' 
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ak  (Ak'kak'  +  Ck'k3k' } 

(1.8) 

bk  =  |I(Bk,kak'  +  Dk ' k^k ' ^ 

and  similar  equations  for  the  creation  operators. 

If  the  particle  operators  satisfy  the  Bose  (Fermi) 
commutation  relations,  then  the  quasi-particle  opera¬ 
tors  again  satisfy  the  Bose  (Fermi)  commutation 
relations.  This  is  guaranteed  by  the  unitarity  of  U. 


NOTE:  In  the  following  theorems, the  upper 

sign  applies  to  Fermions,  the  lower 

F-D 

sign  to  Bosons:  ± 


Theorem  1 

If  I.  [ak'ak' ] ±  kk' 

II.  [bk,b£, ] ±  =  6kk, 

III.  [ak,bk, ] +  =  0 

IV.  [ak,ak,]±  =  0 

V.  [bk,bk, ] +  =  0 


(1.9) 


all  hold,  then  all  of  the  following  are  true: 


I. 

[ak'ak' ] ± 

6kk’ 

II. 

[3k^3k, ] ± 

6kk' 

III. 

[ak/ 3k, ] ± 

0 

B-E 


(1.10) 


. 
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IV .  [  oi^  ,  ,  ]  +  —  0 


V.  [3k'3k'] ±  =  0 


We  shall  prove  only  (1.10)  I  and  V.  Proofs  for  the 
rest  are  similar. 

Proof : 


J, 

I .  [  c^k ,  c^k ,  ] 

=  £  {  (A,  a  +B,  b  )  (A*  .a^.+B*  .b1".) 

qq .  kq  q  kq  qM  k'q'  q'  k'q'  q1' 


±  (A*  ,at,+B*  ,bt,)(A1  a  +B,  b  )} 

k'q'  q'  k'q'  q '  kq  q  kq  q'  B-E 


=  l  (A,  A*  ,6  ,  +B,  B*  .  5  ,  ) 

qq .  kq  k  q  qq  kq  k'q'  qq'; 


+  (AkqBk,q,0+BkqAk'q,0) 
^2^2 

l  (A,  A.*  ,  +  B.  B* ,  ) 

g  kq  k'q  kq  k'q 


kk 


v-  t6k'Bk'h 


=  I  Uc.  a  +D.  b  )  (C.  ,  .a  .+D.  ,  .b  ,) 


qq' 


kqq  kqq  k  '  q  '  q  '  k'q'  q 


±  (C,  ,  ,  a  ,  +D.  ,  ,  b  .  )  (C.  a  +D.  b  )  } 

k'q'  q'  k'q'  q'  kq  q  kq  q 


F-D 

B-E 


=  l  (C,  C,  ,  ,  [a  ,a  ,  ]  .+  D.  D.  ,  ,  [b  ,b  ,  ]  , 

u  ,  kq  k'q'  q  q'±  kq  k'q'  q'  q'  ± 
^2*3 


+  D,  C.  ,  ,[b  ,a  ,  ],+  C,  D,  ,  ,  [a  ,b 

kq  k'q  q'  q'  ±  kq  k'q'  q'  q'  ± 


=  0 


Q.E.D 


* '  pp  p  1  "pp  *'3>  j»i  'p*t  , 


. 


We  shall  also  need  the  following  results: 


Theorem  2 

If  Jf=  l  ^  akak  +  ^  Ak  } 
k  k 

then 

I  •  —  ^k  '  ^k  ' 

.  (2 ) 

II.  3k,^=  ^3k,  +  3k,  . 


(1.11) 


Proof:  We  will  prove  I  only.  The  proof  for 


II  is  the  same. 


ot. 


)(ak'ak)ak  +  Ak  ^kak'V 


F-D 

B-E 


“  I  ^Xv  (  Vv  !  +  a’  a 


_  _  t 


■  Kt 


F-D 


k  v  ~kk 1  T  '~k~k,,~k  ^  “k  pkMk  k,J  b-E 


=  +  ^ak« 


Corollary 

If  is  as  above,  then 


I .  ak , e 


-iH  t/h  -i.34  t/h 
1  o 


ak-  e 


•  ,(2) 


TT  d  -i#t/K  _  -iJ+t/fi  o  e  lXk't/fi 

II.  6.  , e  =  e  e 


(1.12)' 


.. 


' 


Proof:  We  shall  prove  I  only.  The  proof  for 


II  is  the  same. 


ve 


-iJalt/fi  _  V  (-it/fi) 


n 


=  l 


n=0 


n! 


n 


But  ak ,  ak ,  +  Xkl  ^  ak ,  n  1 


n  )  n 

=  )  a 


ak ' e 


-iJt-  t/fi  _  -iJ^t/fi 


=  e 


ak '  e 


-ix/’t/K 


Q.E.D 


To  get  the  time  evolution  of  the  system,  first  calculate 
e.g.  : 


iv^t/fi  -  i  5^  t/fi 

e  a,  e 

k 


=  e 


iJ+t/fi  I  (Ak-kak-  +ck-kek')e_i^t/K 

k ' 


k . 


-i ^  t/fi 

^  ,  (Ak '  ke  ak  '  +  Ck '  ke 


3k. ) 


=  l  [pkk,(t)ak'  +  Qkk,(t)bk,] 

\r  I 


(1.13) 


where 

-iAk^t/fi 

Pkk 1  I  (Av »vAv »v»e  » i- »  e 


. ,  (2)  .  ^ 

-iAk„  t/fi 


k" 


k  "  k  k  "  k 


'k"k  k"k 


) 


-iAk^  t/fi 

Qkk,  (t)=  I  (A* „v  , e  +  C* „,_D 


k" 


k"k  k"k 


k"k  k"k ' 


-iXj/t/fi 
e  ) 


(1.14) 


i  i 

and  similar  expressions  for  k>k,  ak,  and  bk.  This 


enables  us  to  calculate  the  time  dependence  of 


the  number  density  e.g.  of  particles  a  in  state  k: 


n^1^  (t)  =  Tr (a^aRpt)/Tr (pt) 

ti&t/h  .  -iJEit/K 
=  Tr  (e  akak  e  PQ)/Tr(pt) 


Tr  {  * 
qq ' 


x 


[pkq-  (t)aa’+  Qka'  (t)bq,]po}/Tr(po) 


kq 


■  l 

qq' 


Tr  (a^a  ,  e  °) 

P*  (t)P  (t)  _ _ 

l*  kqlt,Fkq'  1  '  -3  if. 


Tr  (e  °) 


Tr(b+b  ,e  3^o) 

+  <q(t)Qkq'(t) 


^WW. 


Tr  (e  ) 


+  Pkq(t)Qkq'(t) 


+  -B^o 

Tr (aTb  , e  ) 

q  q 


-3£f 


o. 


Tr  (e  ) 


+ 

Tr (b  a  , e  ) 

+  Qkq(t)pkq-(t)  - a--er5T - 

Tr  (e  ) 


(1.15) 


-3if_ 

where  p  =  e  is  the  equilibrium  statistical 

o 

operator  for  t  <  0.  The  traces  can  be  easily  evaluated 
using  the  commutation  relations  (1.9)  and  the  corollary 
to  Theorem  2  (replacing^  by  a/3  by  a,b;  X  by 

e(l),  and  it/tf  by  3).  The  result  is 


■ 


9 


t  ~mo 

Tr(a.^ _ L=  6  n^U), 

kk'  1  k  ; 


(1.16) 


Tr (e  °) 


-3  SL 

Tr(b^'e _ L=6  n(o)(e(2), 

-3  kk'  1  k  ' 

Tr  (e  °) 


(1.17) 


t 


Tr(akbk,e  )  =  Tr(bkak,e  )=0  (1.18) 


where 


(o)  ,  (i) . 
n  (£k  > 


Rr  (i) 

-By  P  k  .  . 

e  e  ±  1 


F-D 

B-E 


is  either  the  Bose-Einstein  or  the  Fermi-Dirac  distri¬ 
bution  function,  depending  on  the  statistics  the  oper¬ 
ators  obey. 

Using  (1 . 16 ) - (1 . 18 )  ,  eqn.  (1.15)  for  the  time 
evolution  of  the  number  density  of  the  a  particles 
becomes 


n^1’  (t)  =  l  [|Pkk,  (t)  r  n  t'°>  (^f')  +  |Q  ,  (t)  Iz  n  ^ol  (e^t)] 

k ' 

(1.19) 


2  _(o)  ,„(1) 


2  _  (o)  /e.  (2) 


The  general  structure  of  the  coefficients  iP^^tt) 


i  1 2 

and  |Q,  (t)  |  is  as  follows: 


Pkk '  (t) 


(1) 


P.  .  ,  +  I  P'  ,  cos  (X  v  y  ) 

kk'  L  ,  kk  qq  q  q' 

4/4 


(2),  t 
ft 


+  y  p"  ,  cos (X  ,  -  \KL)) 
i  ,  kk'qq*  n  n 


q^q 


(i)  ,  (i) ,  t 

ft 


r  _  m  i  /  ->  (  2  )  •>  ( 2  )  .  t 

+  )  P i  i  |  ,  cos  (X  ,  -  X  )  zr 

L.  .  kk '  aa '  a '  a  ft 


q^q 


,  kk'qq 


(1.20) 


OI.I) 


M- 

.  (X)  .  (o)  ,  (  s'y6>):lT 

(  „3)  a  lj(j(a  -  — 


.&a- 


9)  iT 


(VI. I) 


,■*«-  t 

,<£)  ,  (o)„  *  _  1  9'Atd):lT 

-J)  n  -x*4 - ?CF - 


9)  IT 


.Ra-  421a- 

(81. 1)  0  =  (°  9,^6^d)^T  =  (°  9 ,  ^d.^6)  ’ll 


9^9riw 


a-'i 


X±  *  >«% 


-  ( (1)  3) 

(f3  n 


' 

-i9qo  9rfd  aoidaidBda  arid  no  pnibn9q9b  ^noidonnd  noxdi/d 

•Y^do  anode 

smxd  add  nod  (5I.I)  ,np9  ,(8I.I)-(dI  I)  pniaU 

39loidiBq  b  9rid  do  Ydian9b  xsdnujn  odd  do  noxdulov9 

8orn.)osd 


[('?  Is)  'o)n  U|(i)fjU0l  +  ((^a)  (o;  i  S|(i),_q|]  J=(:»(I>n 


**v'  '  '  'X~'  “  1  VJ/  '^bT'  J  w/  X 


A\ 

(ei.i) 

1  '  I;:  1  adn9xoi:^d9oo  9fid  do  9nodonnds  .  6n9fl9p  9dT 

_ r  r-i 


d  ,(£) 

*  1  y 


(i) 


pp'^  , _ ^ 
p » p 


d  .  (I)  ,  (I)  „ 

■5L  ;  X  “  ,  X)  aoo  ,  ,  ,  rq  + 

n  P  P  'PPV**  .p^p 


(OS. I) 


d  ,(S) 


■J  /  V  ^  I  /  (  S.  )  . »  in  7 

*  <  PX  "  /)ao°  'PP'^q  i  + 


:  '  ’  » M>rQ ! 

bnj6 

|(i), 

x**! 

Here  P^kl  measures  the  amount  of  a  particles  still 
present  in  the  new  steady  state  as  t  -*  00  .  The  last 
three  terms  give  the  momentum-state  mixing  of  parti¬ 
cles  a  and  b  together,  and  of  particles  a  among 
themselves,  and  of  b  among  themselves,  respectively. 

All  these  expressions  simplify  considerably 
if  we  assume  a  constant  of  field  V^k,  =  Vo6^^,  , 

because  then  no  mixing  of  momentum  states  can  occur. 
In  particular,  we  get 


(1)  ...  (o)  .  (1) .  ,  1 

nk  (t)=n  (£k  l  +  2  J 


V 


o 


2 

v  +  z,f 

o  k 


[  1  -  cos ( X 


(D  _  i  (2)  i  t, 
k  k 


r  (o)  ,  (2)  ,  (o)  ,  (IK  , 

x  [n  (e^  )  -  n  (ek  )  ] 


(1.21) 


with 


C  -  r  -  ^ 


■  -  i(6U>  -  E|2') 

’k  2  k  k  ’ 


The  single  particle  energies  are  given  by 

(D  _  1i2k2 

ek  2m^ 


(1.22) 


where 


*2V2 

=  *Jl_  +  £ 

2m2  o 

the  threshold  energy  of  formation  of  par¬ 


ticles  b. 


01 


Hxds  asIoxdxBq  x  do  dnuoms  arid  a^iusBem  *f^  .*3  aieH 
danl  ariT  .  00  «-  d  ss  9dBda  vbBala  wan  arid  nx  dn9S9:iq 
-xdiBq  do  pnxxxm  9dBds-“mjjdn9fnom  9rfd  avip  amxad  99irid 
pnoms  b  aloxdifiq  do  5hb  *i9/id9pod  d  bns  b  solo 
.  Yl9Vid09q29x  ,  a^v-Xoamsri.!  pnomB  d  do  5xb  ,  aevlsamsrid 
Y  dfiisbxanoo.  u  t>  ;i  »  snoxas9iqx9  939rfd  IIA 

x  {>icV  , :.,4V  '">/  id  do  in.bc  oo  b  9/nussB  9w  dx 
.  xjjooo  nr  o  aadfida  uiubna/nom  do  pnxxxw  on  narfd  ssuBood 

d9p  9w  ,iBluoxdtBq  ni 


V 


[ft  (  !;X  ”  U)ao°  ~  -H  4  +  (  r3)  n  =  (d)  '  X?rn 


(IS. I) 


s — r  s 

A  o 


[((1>a><0,r-  (£b)^n]  X 


rfdxw 


+  oV>S  - 


(I) 


(££.I) 


3  + 


2*~ft  =  (£) 

k 


-XBq  do  noidBimtod  do  bloriassrid  arid  ax  3  an:9/iw 

o 

.d  39loid 
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We  now  specialize  to  a  system  consisting 
initially  at  t  <  0  of  only  particles  a  in  equilibrium 
at  a  temperature  T,  a  constant  external  field  VQ  being 
switched  on  spontaneously  at  t  =  0 .  We  then  get  for 
the  density  of  particles  a  in  the  thermodynamic  limit 


Pi(t>  =  |  I 


(2tt) 


■q  3  i  (1)  ... 
d  k  n£  (t) 


=  p,  (t=0) - ±-r 

1  (2tt)  J 


^d3k  i  n  (o)  (£1|1)  ) 

Vo+5k 


(1) 

lk 


x  [1-  cos (A/"'  -  A/"') 


(2),  t. 
lk 


(1.23) 


Making  the  substitutions  6  = 


y  vo 


m- 


kT  ' 


r  = 


o 

V  ' 

o 


T  = 


2V  t 
o 


y 


m]_  m2 
m2-mi 


and 


x2=  * 


2yV 

o 


(1.24) 


and  assuming  T  is  high  enough  to  use  Maxwell-Boltzmann 
statistics  for  simplicity,  the  time  evolution  equation 
becomes 


Px  (t) 

Px  (°) 


1  -  —  6 

/? 


3/2 


2j  r.  2  0  0  ^1 

x  dx  _-6x  {1_cos  T[1+!,;(x2_r)2]  } 


0 


2  2 
l+\(x  -r) 


(1.25) 


Several  examples  of  this  equation  are  plotted  in 
Fig.  1.1,  on  the  next  page.  A  discussion  of  this 
diagram,  taken  from  Reference  [1] ,  will  be  given  in 
Chapter  IV. 
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(ttS) 


(£S.I) 
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Fig.  1.1.  Typical  examples  of  density  oscillations  in 
the  two  channel  model.  (Taken  from  [1]). 
Upper  graph:  6  =  l.;r  =  5. 

Lower  graph:  dotted  line:  6=  1.;  r  =  10. 

dashed  line:  6=  0.1;  r  =  1. 
solid  line:  6=  1.;  r  =  1. 
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CHAPTER  II 


THE  THREE  CHANNEL  SYSTEM 


§2.1  The  Time  Evolution  Equation 

We  consider  three  kinds  of  particles  a,  b  and 
c  (either  all  bosons  or  all  fermions)  which  can  be 
transformed  into  each  other  via  external  fields. 

Since  the  two  channel  system  was  already  done  in 
detail,  we  shall  proceed  rather  quickly. 

The  hamiltonian  for  this  system  in  creation- 
and  annihilation-operator  representation  is 


l  +42)btK  +  43)c!cJ 

,  KL  K  K  JC  K  K  KL  K  K 


+  0(t)  ^ }  (akbk' +  bkak,} 


kk 


+  Vkk ' ^  (bkck' +  ckbk ’ )  * 


(2.1) 


where 


V(ij)  =  i 
vkk '  V 


dxV(l%)  e 


i (k-k ' )  •  x 


As  before,  if  (x)  =  is  a  constant  inter¬ 
action  then  =  v^i^(S]cki* 
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II  HaTSAHD 


M3T2Y2  J3KKAH0  33>1HT  3HT 


uru  d  ,B  :3.Coj  ftsq  o  3e:;r|J  trib  .  ::ioo  r 

9d  nco  rioxrfw  (anoimsl  lie  10  anoaod  lie  mridxe)  o 
.ablsi^  lBni9  rxs  HIV  !3riio  doss  odni  bsflr/oianBtr:*- 


X  .£§ 


3nob  ybs oils  ;  tsdaye  IsnnBrio  owd  sdd  9oai2 

t&dl&i  b9QDO^q  I  Isria  &v  .Xisisb 


ai  noidBlnsas^qsi;  aodB-i^qo-noidsI iriinns  bfiB 


(  0t’  <E). 


(ljt6.,d  +  ,.(d.  6)  7  v}  {  ii)e  + 


919rfw 


'  C  V  =  (x)  cx  v  li  .exoisd  aA 


'*1 


r  V  narii  ncuios 
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Diagonalization  and  exact  solution 


In  matrix  notation 

(akibklck>! 


where 


1 - 1 

W 

V _ 

VU2) 

v(13)' 

£  = 

v(12)  + 

e(2> 

v(23) 

(2.3) 

v  (13 ) t  v  (23 )t 

E  ^ 

bk 

c. 


(2.2) 


and  E^  and  are  submatrices  with  elements 


( r ^  =  e  (i)  5  and  V  ^ ^  ^ 

k  °kk'  d  Vkk ' 


Jkk'  ”  ^  —  Vvv.  respectively.  The 

hamiltonian  eqn.  (2.2)  can  be  diagonalized  with 


a  unitary  matrix 


B  C 

E  F 

H  J 

/ 


That  is: 


t. o+l  + 


(VBklV  £ 


where 


U  e  U 


r s» 


(2.4) 


(2.5) 


(2.6) 


' 
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/  •  \  /  •  \ 
and  has  elements  A^k ,  = 

quasi-particles  by  a  linear 

annihilation  operators 


^kk'^k^*  We  introduce 


transformation  of  the 


ak  £,(Akk,ak'  +  Bkk '  bk '  +  ^k^k'* 

k ' 


3k  (Dkk,ak'  +  Ekk ' bk '  +  Fkk,ck,) 


Yk  £,(Gkk,ak'  +  Hkk,bk'  +  ^k^k^ 


(2.7) 


with  inverse  transformations 


ak  ^,(Ak'kak'  +  Dk'k3k'  +  Gk'kYk,) 

k ' 


kv  I  '  k  ^k '  Hk '  kYk '  ^ 


'k  ,L  ,  k '  k  k 

k ' 


(2.8) 


ck  “  ^ ,  ^Ck '  kak '  +  Fk'k\'  + 


and  similar  equations  for  the  creation  operators. 
Because  of  the  unitarity  of  U,  the  quasi-particle 
operators  again  satisfy  the  same  commutation  relations 
as  the  particle  operators.  That  is,  if 


[ak'ak,]±  6kk' 
^ak'bk' ^  ±  =  ° 


F-D 

B-E 


(2.9) 


[  r  i  1  +  0 


etc. ,  then 


' 
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[ak'ak,]±  6kk’ 


t  ak '  ^k 1  ^  ±  0 


F-D 

B-E 


(2.10) 


[ak'ak,] ±  =  ° 


etc . 


To  get  the  time  evolution  of  the  system,  use  Theorem  2 
and  its  Corollary  from  Chapter  One  to  first  calculate 

e.  g. 


i&t  /fi  -iJEit/fi 
e  ak  e 


=  e 


i^t/fi  v,* 


l (Ak’kak'  +  Dk’k3k' +  Gk'kYk' )e  ^ 

k' 


-ix'l’t/fi  -iX^t/K 

=  I  (Ak ' ke  ak '  +  Dk ' ke  6k ' 

k ' 


+  Gk'ke 


Yk.) 


-  I  [Pkk'  (t)ak'  + Qkk'  (t)bk'  +  Rkk'  (t)ck' 1  (2-11) 


where 


pkk- (t) 


-  i  A^,  t/'fi 

^  (^k 1  k^k  1  k" ®  ^ 

*  -iX^U/R 

A  _  JS.  \ 

+  Gk'kGk,k"e 


-iA^l t/fi 


(2.12) 


Qkk"  (t)  "  l,  (Ak,kBk,k" 


“iAk^t/fi 

e  +  Dk '  kEk '  k,,e 


-iA^t/fi 


+  Gk'kHk'k"e 


-iAk?} t/fi 


) 


(2.13) 


'**b  =  4  '  •5JO'jtJ0,J 


p.S  k 


' 


(lj(r 


At)..  * 

,  '  * 


.  “  ***° 
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r  *  -iA^t/ft  -iX^t/K 

Rkk"  (t)  “  (Ak'kCk'k"e  +  Dk,kFk'k,,e 


-iA^t/fi 

+  Gk,kJk'k"e  *  ^ 


(2.14) 


t  ,_  +  t 

k 

This  enables  us  to  calculate  the  time  dependence  of 


and  similar  expressions  for  b^,  c^,  a^,  b^  and  c^. 


the  number  density  e.g.  of  particles  a  in  state  k: 

...  Tr  (afa,  p  )  Tr  (ei#t/fi  a+av  p  ) 

_  (1)  /J_\  -  _ k  k  t  _  _  K  K  o 

nk  “  Tr  ( p"  ) 


Tr  (p  ) 


v  r|p  /j_\|2  (o)/  (l)x  in  /,xi  2  (o)  ,  (2)  . 


k' 


kk 


+  |  Rkk,  (t)  |2n(o)  (e^)] 


(2.15) 


-Bit 


where  p  =  e 
o 


o 


is  the  equilibrium  statistical 


operator  for  t  <  0  and  n^°^  (e^^  )  are  the  B-E  or  F-D 
distribution  functions,  depending  on  the  statistics 


the  operators  obey. 


No  momentum  mixing 

All  these  expressions  simplify  considerably  if 
we  consider  a  constant  field  V^-j  ^  6^^,  because 

then  no  momentum  mixing  occurs.  In  particular,  each 
of  the  sub-matrices  of  U,  i.e.  A,B,...J  are  diagonal. 
This  is  proven  in  the  following  three  theorems. 


' 
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Theorem  3 

If  A,  B,  C,  D  are  any  general  block  matrices 
with  A  and  D  square,  then  [2] : 


A  B 
C  D 


=  A 


D  -  CA  1 B 


Proof : 


Theorem  4 


f 

A 

> 

B 

f 

C 

V. 

D 

< 

A 

B 

=  1 

C 

D 

CA 


-1 


0 

I 


A  0 
0  D-CA_1B 


I  A_1B 


0 


A (D-CA  1B)  I  -  1 


Q.E.D 


If  A, B , . . . J  are  square  diagonal  submatrices, 


then 


ABC 
D  E  F 
G  H  J 


A(EJ-FH)  -  B(DJ-FG)  +  C(DH-EG) 


Proof:  Use  Theorem  3  with  the  matrix  partitioned 


as  shown: 


A_j_B  C 

D  !  E  F 

=  |  A  |. 

G  |  H  J 

\ 

E-DA 

H-GA 


B 

B 


F-DA  1C' 
J-GA_1C 


=  |A(EJ-FH)  -  B(DJ-FG)  +  C(DH-EG) 


Q.E.D 


. 


' 
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Theorem  5 

If  e  is  as  given  in  eqn.  (2.3),  then  its 
eigenvalues  are  found  by  solving  the  equation 
1  XI— e:  I  =  0.  In  the  transformation  matrix  U,  each 

I  I  9 

•v 

submatrix  is  diagonal  if  V^-j  ^  =  V  6,  ^ ,  . 

Proof :  Using  Theorem  4, 

0  =  I Al-e I =  |F  (A)  I  (2.16) 

where  F(A)  is  a  diagonal  matrix  with  elements 


F(A)kk'  =  6kk'u3-x2(ek1)  +  £k2,  +  £k3,) 

+  X(£<1)^2)+£<1)e<3)+e<2)£<3)-V<23)2-V^2>2-v<13)2) 

-  (^1)ef>£<3)+2V(12)V(13)V(23)-£‘1)V(23)2-e<2)V(13)2 

-e^3)V(12)  )}.  (2.17) 


Expanding  the  determinant 


0  =  |F(A)  I  =  n  F  ( A ) 


kk 


F(A)kk  =  0  for  a-*--L  k. 


This  third  order  equation  in  A  has  the  solutions 
A  ,  i  =  1,2,3.  The  eigenvalues  are  arranged  in 

K 


the  order  given  in  (2.6)  and  A 

=  '  • 

kk 1  k  kk 


(i) 


has  the  elements 


. 


. 
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The  columns  of  U 


-1 


malized  column  vectors 


rx 


z 

v.  y 


are  formed  from  the  nor- 
where 


etc.,  which  are  obtained  by 


f  > 

X 

(  \ 
0 

y 

= 

0 

z 

0 

k  y 

solving  the  equation 

(2.18) 


(i) 


If  A  =  ,  we  get  the  following  system 


(^i)-e£1,)*n-v(12,yA-v(13,*.-o 


£ 


-  V  ( 12  ’  -  E.'-'Jy,  -V'~'z„  =  0 


L(i)  _  e  (2>  )y  -v(23L 

k  «.  ;yJl  «. 


(2.19) 


-V<13>x„  -V<23>y,+  (X][i)-ej[3))zJl  =  0 


'£ 


for  each  Jl.  There  are  two  cases:  £  ^  k  and  £  =  k. 


£  ^  k  Then  this  set  of  3  equations  (for  a  given  £) 
has  determinant  ^  0  since  and  are  not 

algebraically  related. 


.*.  x£  =  y^  =  z£  =  0  for  all  £  /  k. 

£  =  k  Since  the  above  equations  all  had  det.  0 
and  yet  |Al-e|  =  O,  so  the  three  equations  for  £  =  k 


must  have  det.=  0. 


-ion  edJ  moi5  bamiol  big  “u  to  anmuloo  ariT 


. 

>:to  >v  'rnaloo  besilem 


=  x 


no  siisjjpa  9rfd  pnxvloa  Yd  bsnxB^cfo  9xs  rfoxrtw  , 


(8  C.S) 


m 9daYs  pnxwollol  9rid  dap  9w  ,  .A  *  A  II 


' 

A  '  A3 


0  =  sS((£)3-U^)  +,Y(eS)v- 


(El) 


A  bns  d  \  A  :298fiD  owd  9:ib  9^sriT 


(A  navi.p  b  iod)  anoiifiLips  E  do  d9a  axrld  noriT  M  ^  5. 

■  ■■"■>,,  


ton  3ig  {  L  a  5ns  A  eottia  0  ^  dfiBnxnrxsdsb  asri 


^  A  IIb  io'i 


0  H  .d9b  bsd  IIb  anoxdsupo  9vodB  9rid  sonia 


xol  anoxdBupa  aairid  arid  oa  ,0  =  I  3-IX|  dev  bns 
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..  at  least  one  of  xk'Yk'zp  is  not  equal  to  zero. 

Thus 

Akk '  =  6kk,xk  ' 

Bkk '  =  6kk'yk  ' 

Ckk '  =  6kk,Zk  ' 

all  formed  with  A  =  A^1^ .  The  other  submatrices 
similarly  are  found  to  be  diagonal. 

Q .  E  .  D . 

If  we  furthermore  assume  that  there  are  only  parti¬ 
cles  a  present  at  time  t  =  0,  then  the  time  dependence 
of  the  number  density  of  the  a,  b,  and  c  particles  is 

n^1'  (t)  =  IP^  (t)  |2  n(o)  (e^1] )  i=  1,2,3  (2.20) 


where 


Pkk)  (t)  |2=  1~2a2V2  a~£1)~2A2G2  (l-f2)-2D2G2  a~f3)  (2.21) 

Pkk}  °-2ABDE(1"f1)"2ABGH(1"f2)"2DEGH(1"f3)  (2  •  22 ) 

Pkk) (t) |2=  0~2aCDf (!-f g) -2ACGJ (1“f2) “2DFGJ (!-f 3)  (2.23) 


with 


cos  (A 


(2) 

k 


t 


cos  (A 


(3) 

k 


t 


cos  ( A 


(3) 

k 


t 

fi 


(2.24) 


t269l 

surlT 

•  *** 

•*** 

aeoiilBmduz  ladlo  ariT  .  '  A  =  A  rfdiw  bsm^oi  I  As 

.Isnopsib  9c{  o  :  bn/zc  .:-  :x  vIij^.A  r.raia 
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si  s9loiiisq  o  bas  ,d’,6  Jfi  ft*  Y^i3fl9b  led/nuri  9ri~  3o 


(OS.S) 


£,sa  =  i 


A  »  (O'  ,  S,  ,  t,  (X) , 


--  U>  ( 


Qiariw 


(ss.s)  (  i-x)Hoaas-(  2-i)HoaA£-(  -x) aaaAS-o  =  "|  d) 


xiXiw 


i  (X)  (6.)  M 


3  , (^)  r  <£)  ^ 
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and  the  unitary  matrix's  submatrices  have 

subscripts  k,k.  i.e.:  AEAkk/  etc.  Note  that 


Fkk  (t=°)  I2 


0 


if  i  =  1 


if  i  =  2,3 


as  was  expected,  and  that  by  the  unitarity  of  U, 

1  l^^l2  =  1  • 

i=l 


i.e.:  Conservation  of  total  number  is  guaranteed  by 
the  unitarity  of  U. 

We  get  for  the  density  of  particles  a,  b,  and 
c  in  the  thermodynamic  limit 


Pi  (t) 


b  l  nkiJ  (t) 


(  2  TT ) 


,3,  (i)  .  N 

d  k  n^  (t) 


(2tt) 


d3k|P^)  (t)  |2  n(o)  (e^1’  )  .  (2.25) 


The  single  particle  energies  are  given  by 


2  2 

„  (1)  _  ft  k 
£k  2m^  ' 

(2)  ft2k2 

Ek  =  2TSJ-  +  £2  ' 

(3)  _*V  +  , 

ck  2m3  3 


overt  asoii  Bmdua  a  '  xixd&m  Y^^inu  srfd  bne 


.o  ^  ^  ,  fA  I  A  :.3.X  .  alqilD.  (us 


evs  =  x  51  o 


x  U  do  9f!l  Yd  derid  bne  xb9doaqx9  ax>w  as 


•  '--‘iw'&i  i 


Yci  b99^nsxBup  ax  ladmun  Ifidod  do  rioxl6vx9ano0  :.9.x 


.U  do  YJx^Blxnu  srfd 

a9loidxsq  do  Ydxaneb  arid  aod  dsp  9W 

' 

dim-.I  oi/nsnYbonnsfid  arid  ax  o 


U)  .  7  lu  . 


<«.s>  .(u<3>  :o!n  ;i  U) 


Yd  aavxp  9ib  aaipxsne  sloxdxaq  9lpaxa  t>riT 


,  r>3  +  - -T  =  ,3 


where  z^,  z^,  are  threshold  energies  for  the  formation 
of  b,  c,  particles  respectively.  To  make  the  final 
time  evolution  equation  dimensionless,  write  z  as 


[s(1) 

I 

u 

(i)' 

£  = 

V(12) 

I 

s(2) 

u 

(2) 

E  V(12 

5\/ 

u(1) 

V. 

u(2) 

s 

(3) 

the  submatrices 

have 

the 

elements 

s(1) 

Kk 1 

^2,2 

.  m  k  x 

2mlV12 

kk’ 

q  (i) 

bkk’ 

,*2k2 

2raiV12 

z . 

+  VX2 

■)6kk- 

i  =  2 

V 


(12) 


(2.26) 


(2.27) 


Ikk' 


=  6 


U 


(i)  V 
kk' 


kk' 

(13) 


V 


TI2T  kk 


i  =  1,2 


z  is  diagonalized  by  the  orthogonal  matrix  U,  given 
in  (2.4).  i . e . : 


A  =  U £  O-^V'12’ 


r  £  i 

«v 

\ 

T  -1 

U 

> 

1^12)  J 

u 

fp( 

=  V  P  =  v(12) 


(2) 


(3) 


(2.28) 


Using  the  substitutions 
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(X)0 
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£,I  =x 


(£i)V  (x) 
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(and  dropping  the  subscripts  of  the  submatrices  of  U, 
and  the  P's)  and  using  Maxwell-Boltzmann  statistics 
for  simplicity,  the  time  evolution  equation  (2.25) 
becomes 
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§2.2  Illustration  Examples 


Case  I  Altering  the  Energy  Scale 


The  frequency  of  any  oscillations  in  the  three 

channel  model  is  proportional  to  the  energies  involved 

in  the  model.  To  show  this,  suppose  all  energies 

involved  are  multiplied  by  a  factor  q.  Noting  that 

(12 ) 

all  energies  are  rationalized  in  terms  of  V  ,  we 
must  make  the  following  changes  in  eqns.  (2.26)  to 

(2.30)  :  e  ■*  qg,  (eq.  V  +qV  )  6  q6,  x  x  , 

<v  ~ 

t  qx,  and  U  and  P^,  i  =1,2,3  are  unchanged.  Eqn. 

(2.30)  becomes 
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(2.31) 


Thus  multiplying  the  energies  by  q  causes  the  oscilla¬ 
tions  to  be  q  times  as  fast. 


Equal  Masses  and  Undamped  Oscillations 

If  m^  ~  m2  =  m3 '  then  using  the  definitions  of 


(2.29)  we  get  that 
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where  1^  =  S^j .  The  matrix  of  rationalized  eigen¬ 
values  is  — (Y2y  A  =  p  where 
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where  y^,  y2,  y2  an<^  u  are  constant  in  x.  The  inte¬ 
grals  (2.30)  for  the  particle  densities  are  then 
trivial.  Note  that  6  drops  out  i.e.:  there  is  no 
temperature  dependence.  We  get 


pi  (t> 

p^ToT  = 


r 


(ll 

0 

-  2 

A 

(A 

B 

D 

f  D 
E 

0 

h 

C 

F 

L  / 

(1  -  cos 


A 

G 

+  A 

B 

G 

H 

C 

J 

'd 

fG' 

+  D 

E 

G 

H 

F 

t,  J 

J 

(1 

(1 


cos 


cos  io^t) 


(2.34) 


rIT 


1 

on  8 i.  93arli  :  . 9 •  i  iuo  aqoab  l  dioM 


a 

27 


where 


Thus  the  oscillations  are  undamped.  This  is  not  true 
in  general  because  usually  the  frequencies  are 
functions  of  x,  resulting  in  damping. 


Case  II  Altering  the  V 


(23) 


interaction  energy 


We  set  m-j^  =  m2  =  m3,  r2  =  r3  =  0  and  v^12^=V^13^. 

Thus  the  a,  b  and  c  particles  are  chemically  identical 

i.e.:  the  same  masses  and  no  threshold  energies 

between  them.  is  a  measure  of  the  probability 

of  an  i  into  j  or  a  j  into  i  particle  transformation. 

(12)  (13) 

Since  V  =  V  and  initially  all  particles  are  of 

type  a  (we  shall  also  say  they  are  all  in  "state  1"), 

we  expect  that  p2  (t)  =  p3  (t) . 

From  Case  I,  we  guess  that  the  oscillation 

(23 ) 

frequency  will  increase  as  V  increases.  Physically 

(23 ) 

this  is  reasonable  since  the  larger  V  is,  the  more 

likely  transitions  between  states  2  and  3  are,  and 
hence  the  faster  the  oscillations.  The  frequency 
can  be  found  exactly: 
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Solving  |  UI— M |  =  0  we  get 


2 

2 


/ 


The  eigenvector  for  y^  is  (0,  +  1//2T,  -1//2) . 
Thus  A  =  0  and  the  only  frequency  that  enters  is 


(2.37) 


As  can  be  seen  from  Fig.  2.1,  the  amplitude 

(23 ) 

of  the  oscillations  also  decrease  as  V  increase 


increases 


This  coupling  of  amplitude  and  frequency  is  typical 
for  oscillations  in  non-linear  systems. 

Case  III  Altering  the  threshold  energy 

We  again  assume  that  m]_  =  m2  =  m3  so  t^ie  time 
evolution  of  densities  is  governed  by  eqn.  (2.34). 
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The  eigenvector  for  ] is  {1//2,  -1//2,  0)  . 

Thus  C  =  0  and  for  the  only  frequency  that  enters 
is 

w3  =  — 2 =  y/(r3"1)  +  8  •  (2.39) 

For  r^  >  1,  w3  increases  as  r^  increases.  Assuming 
the  period  and  amplitude  are  coupled  as  before,  we 
expect  that  as  r^  increases  (r^  >  1) ,  the  period  and 
thus  the  amplitude  of  the  oscillations  for 
decrease.  This  trend  is  shown  in  Fig.  2.2.  Phy¬ 
sically,  this  means  that  as  r^  increases,  less 
particles  have  the  threshold  energy  required  to  become 
C  particles  and  so  the  density  oscillations  are 
smaller.  For  r^  < 1;  as  r^  increases,  decreases. 
Fig.  2.3  shows  that  the  amplitude  and  period  are 
coupled  so  that  the  amplitude  of  actually  increases 
as  r^  increases. 

In  any  case,  for  r^  =  0  and  r^  >  1  or  0  <  r^<  1/ 
the  p1  and  p2  density  oscillations  are  greater  than 
the  p3  density  oscillations.  Physically,  this  means 
that  the  p^  oscillations  are  inhibited  by  a  positive 
r^  threshold  energy. 
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Fig.  2.3.  Equal  mass  systems  with  =  0,  0.2  and  1. 


or .  a  ±  ‘:9. 


Case  IV  Long  period  oscillations 


In  the  two-channel  model  the  density  oscilla¬ 
tions  are  due  to  the  factor 

cos[  u(2)-x(1)  )£]  =  cos[/l  +i-(x2-r)2  ^ - -  t]  .  (2.40) 

(12) 

Thus  the  maximum  period  possible  is  T  =  h/2V 

occurring  when  x2  =  r.  For  example  if  V^12^  ^  10  eV 
-15 

then  T  ^ 10  seconds.  However  in  an  n-channel  model, 

as  n  increases  the  periods  can  become  much  longer 

due  to  beats  or  to  closely  lying  eigenvalues.  This 

is  only  possible  for  n  >  3.  The  reason  is  that  in 

(12) 

a  2  channel  system,  Vv  causes  the  smaller  eigen¬ 
value  to  be  less  in  value  than  the  smaller  diagonal 
element  and  the  larger  eigenvalue  to  be  greater  than 
the  larger  diagonal  element.  Due  to  this  widening, 
the  2  channel  system  has  a  maximum  period.  This 
argument  can  be  applied  to  the  n  channel  system. 

By  rotating  two  axes  at  a  time  (method  of  Jacobi 
rotation)  until  the  hamiltonian  is  diagonalized,  we 
see  that  the  smallest  diagonal  element  must  decrease, 
and  the  largest  must  increase.  But  it  is  possible 
in  the  3  channel  system,  for  the  intermediate  eigen¬ 
value  to  lie  close  to  one  of  the  other  two,  causing 
a  long  period  oscillation. 

An  example  is  the  case  where  r 2  =  2,  r^  =  4, 
mi  =  m2  =  m3,  U1  =  1,  U2  =  3.  Then 
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1 


1  2  3  +  x2I  =  M  +  x2I  . 

13  4 


(12)  § 


(2.41) 


V 


Solving  |  yl— M  |  =  0  we  get  \i±  =  0,  y2  3  =  3±  Sl2 .  Thus 
the  three  periods  c*re  T^^.62tt,  T2^  8.62  7t  and 

^  .58  tt  .  It  can  be  shown  that  the  coefficient  of 
the  cos  o>2 t  factor  is  relatively  large.  Thus  the 
density  oscillations  consist  of  a  large  amplitude- 
long  period  oscillation  with  a  small  amplitude-short 
period  oscillation  superimposed  on  it,  as  shown  in 
Fig.  2.4(b).  Figs.  2.5  and  2.6  show  systems  with 
parameters  slightly  different  from  those  of  Fig.  2.4, 
and  with  the  long  period  absent. 

Systems  with  particles  of  unequal  mass 

If  the  masses  of  the  a,  b,  and  c  particles 
are  not  all  equal,  the  time  evolution  equations 
are  much  more  complicated  than  those  for  the  preced¬ 
ing  examples.  The  unitary  matrix  U  and  the  eigen¬ 
values  are  functions  of  the  integration  variable 
x,  so  we  must  use  the  general  equation  (2.30)  rather 
than  eqn .  (2.34). 


It  is  still  possible  to  make  approximations 


to  (2.30)  for  small  times  and  theoretically  at  least, 
for  large  times.  For  t^O,  the  arguments  of  the 
cosines  are  small  and  the  cosine  factors  oscillate 


35 


(a) 

R2  =  1.00 

R3  =  3.00 

M1/M2  =  1.00 
Ml/ M3  =  1.00 
IJl  =  1.00 

1)2  =  3.00 

DELTA  =  1.00 


(b) 

R2  =  2.00 

R3  =  4.00 

M1/M2  =  1.00 
M1/M3  =  1.00 
U1  =  1.00 

U2  s  3.00 
DELTA  =  1.00 


(C) 

R2  =  2.50 

R3  =  4.50 

M1/M2  =  1.00 
M1/M3  =  1.00 
U1  =  1.00 

U2  =  3.00 

DELTA  =  1.00 


Fig.  2.4.  Long  period  oscillations  in  the  3-channel 
system. 


■ 


36 


(a) 

R2  a  2.00 
R3  =  4.00 

M1/M2  =  1.00 
M1/M3  a  1.00 
U1  a  1.00 
U2  a  1.00 
DELTA  a  1.00 


(b) 

R2  a  2.00 
R3  4  4.00 

M1/M2  4  1.00 
H1/N3  4  1.00 
Ul  a  1.00 
U2  a  3.00 
DELTA  4  1.00 


(c) 

R2  =  2.00 

R3  4  4.00 

M1/M2  a  1.00 
M1/M3  4  1.00 
Ul  a  1.00 
U2  a  5.00 
DELTA  a  1  .00 


Fig .  2.5 


37 


(a) 

R2  =  3.0C 

R3  =  4.00 

M1/M2  =  1  .00 
MI/M3  =  1  .00 
ui  =  i.oo 
U2  =  2.00 

DELTA  =  1.0C 


(b) 

R2  =  1 .00 

R3  =  2.00 
M1/M2  =  1.00 
M1/M3  =  1  .00 
Ul  =  1.00 

U2  =  1.00 

DELTA  =  1.00 


(c) 


R2  = 

2.00 

R3  a 

4.00 

M1/H2 

=  1.00 

M1/M3 

=  1  .00 

Ul  = 

O.OC 

U2  = 

3.00 

DELTA 

=  1.00 

Fig.  2.6 
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slowly  since  the  P^ ' s  are  of  order  1.  The  matrix 

elements  A,  B,  ...J  are  also  of  order  1  so  we  can 

assume  that  the  maximum  contribution  to  the  integral 

2  -Sx^ 

comes  from  the  maximum  of  the  factor  x  e  ,  which 
is  at  x  =  1//S".  Thus  the  densities  are  again  described 
by  eqn .  (2.34)  but  with  U  and  the  ' s  evaluated  by 

eqns.  (2.32)  and  (2.33)  at  x  =  1//S”.  Fig.  2.7 
compares  this  approximation  (a)  to  the  exact  solution 
(b)  for  a  typical  example.  Note  that  now  6  (tempera¬ 
ture)  does  not  drop  out  but  determines  where  the 
integral  has  its  maximum  contribution. 

For  t  large  the  cosines  generally  oscillate 
rapidly  over  the  range  of  integration  and  the  above 
argument  fails.  In  this  domain  Kelvin's  Stationary 
Phase  Argument  [3]  states  that  the  maximum  contribution 
to  the  integral  comes  from  the  region  where  uk(x)  = 
p_.-p^/2  has  a  stationary  point  giving  a  1//T  damped 
time  dependence.  The  other  regions'  contributions 
interfere  destructively.  If  there  are  no  stationary 
points,  the  maximum  contribution  is  from  the  endpoints 
of  the  integration  interval  giving  a  1/t  or  faster 
damped  time  dependence. 


(23) 

Case  V  Altering  V  with  unequal  masses 

A  beautiful  example  of  this  small  and  large 
time  evolution  is  a  system  very  similar  to  that  shown 


. 


. 


39 


R2  =  3.00 

R3  =  5.00 

M1/M2  =  0.00 
M1/M3  =  0.00 
U1  =  1.00 

U2  =  3.00 

DELTA  =  1  .00 


TflU/PI 


Fig.  2.7(a)  Approximate  solution  (evaluating  U  and  P^, 

i  =  1,2,3  at  file  point  of  maximum  contri¬ 
bution  x  =  6"'1)  . 


R2  = 

3.00 

R3  = 

5.00 

M1/M2 

=  0.00 

M1/M3 

II 

O 

o 

o 

U1  = 

1  .00 

U2  = 

3.00 

DELTA 

=  1  .00 

Fig.  2.7(b)  Exact  solution. 
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in  Fig.  2.1.  This  time  we  let  m^/ir^  =  m^/m^  =  0, 

=  6  =  =  1  and  we  again  vary  U^.  This 

system  is  similar  to  that  in  Fig.  2.1  because  at 
the  point  of  maximum  contribution  (x  =  1/ /&)  their 
hamiltonian  matrices  (2.36)  are  the  same.  However 
this  time  we  expect  damped  oscillations  because  of 
the  unequal  masses. 

The  oscillation  frequency  can  be  found 
explicitly.  The  hamiltonian  matrix  is 


(  £ 

X 

i— H 

Ul] 

x  -1 

1 

i — 1 

1 — 1 

r2 

U2 

= 

1 

0 

U2 

1 — 1 
D 

U2 

r3> 

1 

U2 

0 

M  +1 


(2.42) 


The  solutions  of  the  equation  |XI  — M  |  =  0  are 

2 

X  *“  (  1—U  ~  )  -1  ———————————— 

A1  =  ~U2'  A2  3  = - 2 -  ±  j/(x  -(1+U2)  )  +  8  .  (2.43) 

The  eigenvector  corresponding  to  is 
(0,  +1//2,  -1//Z)  so  A  =  0. 

Thus  the  only  frequency  that  enters  is 

u3  =  |  /(x2-(1+U2))2  +  8.  .  (2.44) 

2 

(jo ^  has  a  stationary  point  at  x  =  1  +  so  the  period 

for  asymptotically  large  times  is  T  =  2 rr/oa  =  /2  ir.  For 
small  times  the  frequency  is 


, 


03 


3 


x=l/  /& 


8 


(for  6  =  1)  and  the  period  is  T  =  4 tt/  v^u ^  +  8.  Thus 
we  expect  initially  a  period  dependent  on  U2  (longer 
for  small  U2 ,  and  shorter  for  large  U^)  and  finally 
for  large  t  a  period  of  /2"  tt  independent  of  U2 . 
Furthermore  the  oscillations  should  be  damped  as 
l//r.  This  behaviour  is  illustrated  in  Fig.  2.8. 

In  2.8(a),  T  =  /2  tt  throughout.  In  2.8(b)  the  period 
increases  from  T  =  -j  tt  to  T  =  /z  it.  And  in  2.8(c) 
the  period  dramatically  increases  from  T  =  4 tt/  /24  s 
.  8  tt  to  T  =  it  at  about  time  t  =  2  tt  or  3  tt  . 


Case  VI  Altering  the  c-particle  mass,  m3 

As  a  final  example  we  study  the  trend  in  the 
density  oscillations  as  the  c  particle  mass  changes. 
We  choose  the  parameters  U-^  =  U2  =  r2  =  r^  =  6  =  1 
and  m1/m2  =  0.  Since  eqn.  (2.30)  in  this  case  cannot 
be  solved  analytically,  we  again  make  the  approxima¬ 
tion  that  the  maximum  contribution  to  the  integral  in 

-1/2 

(2.30)  arises  from  the  point  x  =  6  /  .  Then  the 

hamiltonian  equation  is 


M  +  I 


(2.45) 


where 


' 


- 

:  t-t  ■  U  »aooMo  oW 


; 
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g).00  1.00  2.  CO  3. CO  4.00  S.CO  6.00  7.00 


Fig.  2 . 8  (b) 


1 — 1 
e 

m 

m2 

m 

=  0.;  r2  =  r3  =  6  =  Ux  =  1 .  ;  U2=  1 


0.00  0.20  0.40  0.50  0.80  1.  0.00  0.20  0.40  0.60  0.80 
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SJ.OO  1.00  2.00  3.00  H.00  5.00  6.00  7.00 


O 

CD 


O 

CD 


CD 

3* 


o 

c\j 


O 

CD 


• — l 

s 

_  mi 

m2 

m3 

Fig.  2.8(c)  —  =  —  =  0.;  r2  =  r3  =  6  =  U,  -  1 . ;  U0=  4. 


§5.00  1.00  2.00  3.00  4.00  5.00  6.00  7.00 


M 


1 

1 


(2.46) 
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0 

1 

1 

v 


1 

0 

1 


m 


and  m  =  .  The  solutions  of  the  equation  |  y I  — M  |  =  0 

are 


vi =  _1 

H2  3  =  ^  ±7  J (m_1)2  +  8  •  (2.47) 


The  eigenvector  corresponding  to  y^  is 
Thus  C  =  0  and  the  only  frequency  that 


3 


(1//2,  -1  //l,  0)  . 

enters  for  p3  is 
(2.48) 


with  the  corresponding  period 


(2.49) 


Due  to  the  stationary  point  at  m  =  1  the  asymptotic 
period  for  p3  is  T  =  /2  tt  regardless  of  initial  condi¬ 
tions.  If  we  again  guess  that  the  period  and  amplitude 
are  coupled,  then  we  expect  the  period  and  amplitude 
of  the  p3  oscillations  to  increase  as  m  increases  from 
0  to  1,  and  to  decrease  as  m  increases  for  m>  1.  This 
trend  is  borne  out  in  Figs .  2 . 9  (a)  to  (d) . 

Another  interesting  point  is  that  in  the  limit 


of  large  m  the  eigenvalues  are  y1  -~1,  y2  -  If  and 
y 3  =  m.  The  eigenvector  corresponding  to  y3 


m. 


in  this 


' 

iio<3  erf*  o3 
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TflU/PI 


Fig.  2.9(a) 


m2  m3 


0 


1. 
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TRU/PI 


Fig.  2.9(c) 


r2  r3  !•? 


m. 


m. 


=  0.; 


TflU/PJ 


Fig.  2 . 9  (d)  r3  =  r3 


0 
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limit  is  (0,0,1)  and  thus  G  =  H  =  0.  So  for  and  p2 

the  only  frequency  that  enters  for  large  m  is  = 
(y2”V1i)/2  =  1.  For  p^  the  limiting  frequency  is  of 
course  =  (y3-y2)/2  =  (m-l)/2.  These  trends  are 
illustrated  in  Fig.  2.9(d). 


§2.3  The  n-Channel  System 


We  consider  n  kinds  of  particles 
av  (all  bosons  or  all  fermions)  which  can  be  trans¬ 
formed  into  each  other  via  external  fields. 

The  hamiltonian  for  this  sytem  is 


n 


= 


l  l  J^a^a^ 

i=l  k  k  k  k 


+  6(t).l  J,  Vkk'J  (aki>+4')  + 

i<]  kk' 

+  a]^j)+a^t))  •  (2.50) 


If  we  assume  that  the  external  fields  V 

(ij)  _  tt  ( i  j  ) 


(ij) 


constant,  then  V^,  =  V 

i  <  j .  Then 


are 

6kk.  for  i,j  =1,2,.. 


•  n; 


M  =  l  {  l  ek1>akl)+ak1>  +  e(t)  l  V,lj)(aW^'  + 

k  x=l  i<j 


n 


I  (  l  a^(l)  (E,J^aij)  } 
k  i,j=l 


k  ij  k 


(2.51) 


where 


(Ek)ij 


_  „  (i) 


=  £ 


1 J 


+  (1- 


6ij)V 


(ij) 


(2.52) 


This  hamiltonian  equation  can  be  diagonalized  with 
a  unitary  matrix,  U.  That  is 

u  =  l  {£  “ki)+“ki)*ki)}  <2'53) 

k  i 


where 


■ 


. 
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*-(l)  =  I  (U,  )  .  .  (E,  )  .  (U* )  . 


k 


k  ij  k  j  £  k  i£ 


(2.54) 


and  where  we  have  introduced  quasiparticles  by  a 
linear  transformation 


a 


(i) 


=  l  (Ik)  .  .aj^ 

h  kink 
D 


(2.55) 


with  inverse  transformation 


aki)  - 


(£) 


(2.56) 


The  quasiparticles  obey  the  same  commutation  relations 
as  the  particles.  That  is,  with  no  momentum  mixing, 
if 


[a 


a<3>\  =  6.. 

k  '  k  ±  l  j 


F-D 

B-E 


(2.57) 


[a'1’,  aW]±  =  0 


then 


r  (i)  (j)t-.  _  * 

[oik  -  ak  ]±  -  hj 


F-D 

B-E 


(2.58) 


r  (i)  ( j )  i  _  n 

[ak  ,  ak  ]±  -  0  • 


For  example, 


(i)  _  ( j  )  t 


akJ  1 1 1  =  l  <ViW  (uk>jA 

£m 


(£)  / tt*n  „  +  (m) 

j: 


.+  (m) 


±  l  <Vjmak .  (Uk>Uak 


U) 


£m 


. 


=  Jm(VuK>jmIak1)'ak(m)l± 


6.  . 
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Theorem  One  and  its  corollary  from  Chapter  One  are 
unchanged  for  n  particles.  That  is 


ak1)e" 


W1’  +Aki)“ki> 


i^t/«=  (i) 

LX 


(2.59) 


-iX^’t/K 


Thus 


(2.60) 


ijcit/h  (i)  -ijlt/h_ 


=  l  (Uk}jie 

3  J 


an 


=  IK’ji e 
3 


-iXj^V* 


-iA  [)  ^  t/K  ,  0  v 

=  Yl  K’ji'Vj*  6  k  >a‘*> 

^  3 


=  ItP^af  > 


(2.61) 


This  enables  us  to  calculate  the  time  dependence  of 
the  number  density  of  particles  a^  in  state  k: 


n(i)(t)  _  Tr(aki)taki)pt) 

k  Tr  (Pt)  * 

Tr{  I  <Ua“>  +  Pklm 
im 


Tr  (e 


-e#. 


(2.62) 


0is.  uoL  60  oi  8 
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-mQ 

where  e  is  the  equilibrium  statistical  operator 

for  time  t  <  0.  Using  the  fact  that 


„,+(£)  (m)  Wo, 

Tr  (ak  a^  e  ) 


=  6 


Tr (e  °) 


£m  -3y 


e  e 


±  1 


F-D 

B-E 


(2.63) 


eqn.  (2.62)  becomes 


(i)  ,,  .  ,  v  ^*i£  o  (o)  .(£) 


K  (t)  =  l  p"~  p 

Jim 


k  'w/  ~  Ln  "  k  *  k  ^£m  n  ^£k 


=  IIP 


i  A  .  2  (o)  ,  (£) 


£ 


k 


n 


<£k  ') 


(2.64) 


where 

i„i£  i  2 


|T  (U,)  .(U*)..(U*)  «(U.).n  e 

1  h  k  mi  k  ]i  k  m£  k  ;j  £ 


3m 


=  5ikJlj  +  Pik£jm  cos(Akm)-Xk3>)  I  •  (2-65) 


Finally,  the  time  evolution  equation  for  the  i 


.  th 


species  is 


Pi(t)  =  -^3 
1  (2 it)  J 


d3k  (t) 


(  2  tt) 


f-3.  v  lt5i £  1 2  (o)  /_.(£)  \ 

d  k  >.  Ipv  I  n  (eic  }  * 


(2.66) 


£ 


If  there  are  only  a^  particles  present  at  time  t=0 


then 


. 


' 

(*3.S)  .  .;<  k-rpi ♦  ti  j  - 
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Pi  (t) 


d3k  Ip 


ij 

k 


)  . 


An  example  of  this  type  of  system  is 


Here , 


a(1) 

t 

a<2> 

^  N 

a(3)  ^  a<4> 

e  would  be 

E ( 1 )  V(12)  0 

*\ 

0 

E  (2>  V(23) 

v(24) 

e<3> 

v<34> 

V 

E  ^  ' 

(2.67) 


CHAPTER  III 


ZERO  TEMPERATURE-TWO  CHANNEL  FERMION  SYSTEM 


§3.1  Time  Evolution  Equation 


In  this  chapter  we  derive  the  zero  temperature 
limit  for  a  two-channel  system  with  both  a  and  b 
particles  present  at  time  t  =  0.  We  start  with  equa¬ 
tion  (1.21) : 


„  ( 1 )  _ 

nk  (t)  - 


(o )  ,  (1).  ,  1 

n  (£k  )+2 


V' 


2  2 
V  +£“ 
o  ^k 


[1  -  cos  (A]^1)  -X^2)  )  x 


X  (n^ufW0^)} 


(3.1) 


where 


(o)  ,  (IK  r  3yl 
n  (e^  )  =  [e 


3 


*2,  2 
h  k 

2m, 


+  1] 


-1 


(3.2) 


✓  2,  2 
3(B_L_+  e  ) 

(o)  ,  (2) .  By2  2m2  °  ,  i i  1  • 

n  (e^  )  =  [e  e  +1] 

y  and  u9  are  the  Fermi  energies  of  the  a  and  b  parti- 
cles  respectively.  In  the  zero  temperature  limit  the 
distribution  functions  (3.2)  become 


2  2 

lim  n  (o)  ( e 1  ’  )  =  1-  0  (fjt--  Uj_)  =  1-  6  (k - ^ 


/2mn  y 


1^1 


3"^ 00 


) 


lim  n(o)  (e<2)  )  =  1-  0  ( 

3^o° 


j/2,  2  /2m0  (y  0-e  ) 

fh  k  \  t  2^2  o  x 

+  e  -  y 0 )  =1-0  (k - s - ) 


2nu  ’  °o  K2 


h 


(3.3) 
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where 


'o 

if 

x  <  0 

e  (x)  = 

• 

1 

if 

X 

IV 

o 

Before  proceeding 

let  us 

first  calculate  p,  (0) ,  p9(0) 

and  the  ratio  p2(0)/p^(0). 
1 


Pt  (0)  = 


(2tt) 


3,  (o)  .  (1)  . 

d  k  n  (e^  ) 


(2tt) 


h 


0 


k2dk 


(2m1y1) 

2  T 

6tt  fi 


3/2 


Po(0)  = 


(  2  7T  ) 


-.3,  (o)  ,  (2)  . 

d  k  n  (e^  ) 


(2m,<u.,-e  )  )  3/2 


(  2  TT ) 


,2,,  v  2'K2  o 

k  dk  = 


0 


2  3 
6tt  h 


(3.4) 


p2(0) 

P3_  (0) 


m2  ^2  _  £o  13/2 
—  (— - -) 

ml  U1 


Using  the  abbreviation 

V2 


F(k't>  =  \  2  °~2  [1 ' cos(Xk1)_xk2> >  r] 


Vo  +5k 


(3.5) 


we 


get  from  eqn.  (3.1)  for  the  density  of  particles 


a  in  the  thermodynamic  limit 


■ 
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/2m2 ( y 2  —  e0 ) /& 


Px  (t)  =  Px  (0)  + 


4iTk  dk  F  (k ,  t ) 


0 


/2m1]i1/fi 


4Trk  dk  F(k,t) 


0 


Y. 


Px  CO)  - 


4nk  dk  F  (k , t )  . 


/2m~  (iu-£  T 

2^2  o 


fT 


(3.6) 


Defining  the  parameters  x,  r,  t  and  y  as  in  Chapter 
One ,  i . e . : 


x  = 


lik 


/2yV 

o 


r  = 


o 


V 


o 


T  = 


2V  t 
o 

h 


y  = 


mlm2 

m2_ml 


(3.7) 


the  time  evolution  equation  (3.6)  for  the  a  particle 
density  can  be  written  as 


P1  ^ 


=  P1(0) 


i  - 1  <—  —  > 

2  \i± 


mn  y,  1/2 

(—  — ) 

V  3/2  U  v  ' 


X 


( 


m2  U  eo  1//2 


y  v 


-) 


o 


X 


2  , 

x  dx 


,  1,  2  ,2 

1+  ^-(x  -r) 


[1  -  cos  T 


A  ^  1'  2  72, 

/  1 +  x(x  -r)  ] 


or  finally  as 


'  I 

« 


/ 


- 

'•i  V 


« 
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px(t) 


P1(0) 


Px  (0)+p2  (0)  p'1  (0)+p2  (0) 


x-i  fi3/2 


1//3" 


,  p?(0)  1/3 

7f 


2, 

x  dx 

I  2  7 

1+  ^-(x^-r)Z 


[  1  -  COS  T 


/TTi 


(x  -r )  ] 


(3.8) 


where 


V 

6  =  JL 

mi 


§3.2  Kelvin's  Stationary  Phase  Argument 


To  study  the  eqn.  (3.8)  just  derived  above, 
it  is  convenient  at  this  time  to  discuss  the  evalua¬ 
tion  of  integrals  like  the  one  in  (3.8)  for  asympto¬ 
tically  large  times. 

As  we  have  already  said  in  Chapter  Two,  for 
an  integral  of  the  type 


f  (t) 


3 

► 

g ( x ) co  s [ h ( x ) t ] dx 


(3.9) 


a 

where  g  is  continuous  and  h  is  twice  continuously 
differentiable  we  can  make  approximations  to  it  for 
small  t  and  large  t.  For  small  times  the  argument 
of  the  cosine  is  approximately  constant  over  the  range 
of  integration  and  consequently  the  maximum  contribu¬ 
tion  to  the  integral  comes  from  the  region  x  ^  n 
where  g  has  its  maximum  value.  Thus  for  t  ^  0 

3 

e 


f  (t)  ^  cos  [h  (n)  t] 


g(x)dx 


(3.10) 


a 

For  large  times  the  major  contribution  to  the  integral 
comes  from  the  neighborhood  of  the  endpoints  a  and  3, 
of  the  interval  and  any  points  at  which  h' (x)  =  0; 
the  latter  being  more  important  in  the  first  approxima 
tion  [3],  Contributions  from  other  regions  cancel  by 
destructive  interference. 


. 


' 


. 
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Suppose  first  that  A  is  the  only  stationary 
point,  i.e.:  h'  (A)  =  0,  a  <  A  <  $  and  assume  that 
h" (X)  >  0.  Expanding  h  about  x  =  X 

h(x)  -  h(A)  =  u2  .  (3.11) 


With  this  change  of  variable 

X+e 


f  (t)  ^ 


g  (x)  cos  [h  (x)  t]  dx 


X-e 


(3.9)  becomes 


u. 


-U- 


2u 


g  (x) 

irnin" 


cos [t (h  ( A)  +  u  ) ] du 


(3.12) 


where 


U;L  =  [h(A-E)  -  h  (X)  ) 1/2 


u2  =  [h(A+e)  -h(A)  ] 1/2 


(3.13) 


We  have  assumed  that  only  the  region  where  u  ^  0 
contributes  to  the  integral.  For  u  ^  0  we  may  let 


g  (x)  -*■  g  (X) 


2u 

h'  (x) 


2  1/2 

th"7AT'] 


so  that 


f(t)  ^  lhirTxr1 


1/2 


u. 


g  (X) 


cos  [  t  (u  +  h  (A))]  du 


-u 


(3.14) 


(3.15) 


' 
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By  the  same  argument  as  above  we  may  let  u^  and 
go  to  infinity  and  we  get  after  integration 


F(t)  %  [ 


2  7T 


th"  (A) 


i/2  ^ 

■]  g  ( A)  cos  [th  (A)  +  j] 


(3.16) 


This  is  the  result  of  Kelvin's  Stationary  Phase 
Argument. 

The  contributions  arising  from  the  endpoints 
are  much  more  difficult  to  obtain.  We  will  merely 
state  the  result  given  in  [3] .  Assume  that 

(1)  we  have  broken  up  the  integration  interval 
and  made  changes  of  variable  so  that  h(x)  is  mono- 
tonically  increasing  in  each  subinterval  a  <  x  £  3, 

(2)  g(x)  is  continuously  differentiable  for  a<  x <  3, 

(3)  h(x)  is  differentiable, 

(4)  h' (x)  =  (x-a) P  1(3-x)a  1h1 (x) 

where  p  >  1,  a  >  1  and  h^Cx)  is  positive  and  contin¬ 
uously  differentiable  for  a  <  x  <  3,  and 

(5)  o<  v  <  1,  0  <  y  <  1  where  1-v  and  1-y 
are  the  orders  of  possible  singularities  of  the 
integrand  at  a  and  3  respectively. 

Then  to  first  order 


3 

g 

, 

a 

% 


(x) (x-a) V  1 


k  ( 0 )  p /^\ 
P 


MO) 


(B-x)U  1cos [h (x) t] dx 
cos  [h  (a)  t  +  j  ^-1 

^7p 

cos[h(@)t  -  \  £] 


'V 


a 


(3.17) 


* 


' 


. 


where 


k (u)  =  g1 (x (u) ) u1  V  ^ 
up  =  h(x)  -  h(a) 


(3.18) 


in  the  first  term  and 


£(v)  =  g1(x{v))v1  y 
v°  =  h  ( 3)  “  h  (x) 


(3.19) 


in  the  second  term  and 


g-^  (x)  =g(x)  (x-a)V  1(8-x)'J  1  .  (3.20) 

Notice  that  if  h' (x)  f  0  anywhere  in  a  <  x  <  3  then 
by  assumption  (4)  p  =  a  =  1.  Furthermore  if  g^(x) 
is  bounded  and  continuously  differentiable  in  a  <  x  < 
then  v  =  y  =  1  and  the  contribution  of  the  endpoints 
to  the  integral  is  0(t_1).  On  the  other  hand  the 
contribution  of  any  stationary  points  was  found  in 
eqn.  (3.16)  to  be  ^t 


' 


_ 


§3.3  Example 


With  the  results  of  Section  3.2  in  hand  we 
now  look  at  some  examples  of  eqn.  (3.8).  We  set 
r  =  4,  so  that  (using  the  above  notation) 


h  (x) 


[1  +  i(x2  -  4)2 


1/2 


(3.21) 


h(x)  has  a  stationary  point  at  X  =  2  as  shown  in  the 
right-hand  diagrams  of  Fig.  3.1.  To  illustrate  the 
results  of  Section  3.2  we  take  three  cases  of  eqn. (3.8) 
with  different  integration  intervals. 


Case 

(a)  : 

II 

tO 

< 

a=2 . 5 

< 

3=5 

Case 

(b)  : 

a=l 

< 

>• 

!l 

to 

< 

3=3 

Case 

(c)  : 

1 

a=2 

< 

3=1 

< 

>-> 

ii 

to 

as  shown  in  Fig.  3.1  (a),  (b)  ,  and  (c)  respectively. 

Case  (b)  does  and  cases  (a)  and  (c)  do  not  include 
the  stationary  point  in  their  intervals. 

Note  that  6  must  be  large  to  have  large 
oscillations.  Therefore,  in  light  of  the  fact  that 
we  want  3  =  5  in  case  (a) ,  let  us  interchange  the 
limits  and  change  the  sign  of  the  integral  in  eqn. (3.8) 
Since  case  (a)  is  similar  to  case  (c) ,  we  will  discuss 
only  cases  (b)  and  (c)  at  length. 


■ 
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(b) 

ROl/f  ROl  +-R02  )  r  0.036 
R  =  4  .00 


DELTA  =  1.00 


Fig . 


zero 


3.1.  Examples  of 


temperature  oscillations. 
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Case  (b)  :  To  get  a  =  1  we  let  6=1  and  thus 
3=3  if  p2  (0) /p-^  (0)  =  27.  Applying  result  (3.16) 
with  h(x)  as  above  and 


g  (x) 


2 

x 


(3.22) 


we  get  the  following  results 


(1) 

gU)  =  4 

(2) 

h (A)  =  1 

(3) 

tr 

>• 

ii 

(4) 

3 

• 

g (x) cos [h (x) t ] dx 

^  cos  (T  + 

7T 

4 


P±  (t) 


P1 (0) +p2  (0)  ^  p1 (0) +p2  (0) 


Px  (°°) 


-0.27 


/  ,  TT  \ 

COS  (T  +  j) 

/t 


(3.23) 


This  expression  is  a  very  good  approximation  to 
Fig.  3.1(b)  even  as  early  as  the  second  oscillation, 
verifying  that  the  asymptotic  frequency  comes  from 
the  stationary  point. 

Case  (c)  :  To  have  a  =  ^  and  3=1/  set  6=4 
and  p2(0)/p1(0)  =8.  The  stationary  point  X = 2  is 
outside  the  integration  interval.  Consequently  the 
asymptotic  behavior  of  the  time  evolution  equation 
(3.8)  arises  mainly  from  the  integral’s  endpoints  and 
is  found  using  result  (3.17) .  Following  the  steps  of 


■ 


■* 


the  argument : 


(1) 


x 


g  (x)  =  g1  (x)  = - j - ^ 


1+  J(x“  -  4) 


~2  is  continuously 


differentiable  for  a  <  x  <  3 


g(a=^)  = 


16 

289 


g  ( 3=1)  =  yj 


(2) 


12  2  i/2 

h ( x )  =  -[1  +^(xz  -  4)Z] 


h  (a)  =  - 


17 

8 


h(3)  =  - 


/T3 


(3) 


h'  (x)  =  h1  (x)  =  - 


x(x  -  4) 


l]  1  +  j(x2  -  4)2 


is  positive 


and  continuously  differentiable  for  a  ^  x  £  3 


h' (a)  = 


1_5 

34 


h'  (3)  = 


/I3 


(4) 


p  =  a  =  v=  ]i=l 


(5) 


g (x) cos [h (x) t] dx 


^  {k(0)cos[h(a)x  +  ^]  -  £  (0)  cos  [h  ( 3)  t.  -  j]  }/t 


Tr¬ 


ee) 


To  determine  k(0)  use  eqn.  (3.18) 


k  (0)  =  (a)/^ 


—  9  ( Qj1 )  _  3  2  ^  c  c 

"  h'  (a)  733 


x=a 


(7) 


Determine  MO)  using  eqn.  (3.19): 


MO)  =  g^B)/^ 


x=B 


g(B) 

”  h'  (8) 


3/T3 


^  -.37 


(8) 


Thus 


r  32  r  17  TT, 

^255  COS  8  2  + 


3  /13  2 


r/13  ,  TT,  1  / 

COS  [— T  +  T]  >/t 


The  time  evolution  equation  becomes 


P-l  (t)  P-,  (°°) 

- TnYT - TTvT  ^ - TrfTT - T7rT--{0.167  sin(2.125l) 

P1(0)+P2(0)  P1(0)+P2(0) 


-  0.49  sin (1 . 8 t) }/t .  (3.24) 


The  frequencies  arising  from  the  two  endpoints  cause 
beats  of  frequency 


CO,  -  C09 

CO  =  1  x-  ^  0.1625 


(3.25) 


or  period  T  = 


2tt 

co 


^  12.3  tt.  The  phase  or  average  fre¬ 


quency  is  co  =  oo,+go2/2  ^  1.96  with  corresponding  period 
T^1.02  tt  .  This  value  for  the  phase  period  is  only 
approximate  since  the  amplitudes  in  (3.24)  are  not 
equal  as  was  assumed  to  derive  (3.25).  In  fact  the 
period  should  be  greater  since  the  longer  period 
sinusoidal  has  the  greater  amplitude. 


. 


. 


> 


?  u  Cbv 


All  these  conclusions  are  borne  out  in 


Fig.  3.1(c)  and  Fig.  3.2.  Fig.  3.2  best  shows 
the  beat  phenomenon  and  also  the  fact  that  (3.24) 
is  a  good  approximation  almost  from  time  t  =  0. 
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Fig.  3.2.  Plot  of  Fig.  3.1(c)  for  0  <  x  <  30  tt  (p^COJ/p-^CO)  =  8.  ;  r=4.;6=  4. 

This  shows  the  beats  of  period  T  =  12.3  tt  arising  from  interference 
of  the  endpoint  contribution  to  the  integral  in  eqn.  (3.8). 


CHAPTER  IV 


TWO  TEMPERATURE-TWO  CHANNEL  SYSTEM 


§4.1  Time  Evolution  Equation 


In  this  chapter  we  study  a  two-channel  system 
initially  with  the  a  particle  in  equilibrium  at  tem¬ 
perature  T  =  separate  from  the  b  particles  in 

equilibrium  at  temperature  T^  =  l/kg32  (kg  is 
Boltzmann's  constant).  At  time  t  =  0  the  particles 
mix  and  an  external  field  is  switched  on  spontaneously. 
The  time  evolution  of  the  number  density  is  described 
by  eqn.  (1.21) 


„  (1) 
nk  (t) 


V 


o 


2  ^2 
V  +£, 
o  k 


[1  -  cos (e^1^ x 


(4.1) 


where 


-3, 


n2k2 


,  ,  / >  3,  e  1  2m. 

n(o)  =  e  1  C1  e  1 

B,e  +e  ) 

(o)  ,  (2) .  2  c2  2  2m2 

n v  '  (e^  )  =  e  e 


(4.2) 


We  have  assumed  that  T^  and  Tb  are  high  enough  to  use 

Maxwell-Boltzmann  statistics.  e  and  e  are  the 

1  2 

chemical  potentials  of  the  a  and  b  particles  respec¬ 
tively  and  eQ  is  the  threshold  energy  for  the  formation 


,  .  I  v.' . 


' 
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of  particles  b. 

It  is  convenient  to  first  calculate  p^(0), 
P2(0),  and  p2(0)/p1(0),  where 


P±(0)  = 


(2  tt) 


4  tt  k2dk  n^°^  ( ^  ) 


i  =  1,2 


Using  the  fact  that 


oo 


0 


,  2  -6k 
k  e  dk 


/tt  1 


-r  777 


we  get 


P1(0) 


p2  (0) 


P2  (0) 

pYToT 


3  *i  £  0/0 

1  c-,  m,  3/2 

e  1  ( - i-y) 

2tt3-l^ 

3~  (e  -e  )  o/o 

2  c2  o  m2  3/2 

e  ( - 7) 

2tt32^ 

^2  ^ £c  "£o) 

e  2  V/2 

3l£c,  ml  32 


(4.3) 


(4.4) 


(4.5) 


Using  expressions  (4.5)  calculate 


n<°)(£(2>)  _  n(o)(e(l)) 

k  *2  ,2 

30  (e  -E  )  -3o  4- — 

2  c2  o  2  2m2 

=  e  e  -e 


6l£c 


-B 


*2k2 


1  2m. 


27r3-.fi2  3/2 
=  (  inT  *  Pl(0) 


o  fi2k2 

p2  ^  3/2  ^  1  2ml 

Ip^w  Y  e 


-3 


*2k2 

1  2m, 


(4.6) 
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where 


y  = 


^1  ^2 

m2  B1 


(4.7) 


Thus  the  time  evolution  of  the  a  particle  density  is 
given  by 


Px(t)  = 


(  2  7T  ) 


d3k  n^X)  (t) 
k 


-  P-^  (0 )  + 


(2tt) 


k2dk 


V' 


o 


2  2 

v  +sr 

o  k 


[1  -  cos (A 


X  {n‘°>(£<2>)  -  n<°>  («=“>>} 


(4.8) 


Define  substitutions  similar  to  (1.24)  of  Chapter  One 
That  is 


6  =  —  $,V 
m^  1  o 


e 

_  o 


V 


T  = 


2V  t 
o 


o 


y  = 


mlm2 

m2~ml 


o  ^2,  2 

2  f  k 

x  = 


2  yV 


o 


(4.9) 


Using  the  definitions  (4.9)  and  dividing  by  p1  (0) +P2  (0)  ,• 
eqn .  (4.8)  for  the  a  particle  density  becomes  finally: 


P1(t) 

Px (0)+P2  (0) 


Pl(0) 


P± (0)+P2 (0) 


1  +  6 

/it 


2  .3/2 


2  , 

x  dx 


0 


1,1/2  )2 

1+  ^-(x  -r) 


x 


x  [ 1  -  COS  T 


7l+  hx2-r>'2]  x 


r  3/2  p2(0)  -Y«x2  -6x2 

x{^  ^T07  e  -  e  } 


(4.10) 


I 
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§4.2  Examples 

Before  discussing  any  specific  examples  let 

us  obtain  the  general  behavior  of  the  time  evolution 

equation  (4.10)  for  large  and  small  times. 

For  large  times  we  can  apply  the  arguments  of 

Section  3.2.  The  endpoint  contributions  to  the  inte- 

-3/2 

gral  of  eqn.  (4.10)  are  0(t  '  )  since  the  integrand 

vanishes  at  x  =  0  and  x  =  °°.  Thus  to  lowest  order  we 

need  only  apply  the  Stationary  Phase  Argument  to  the 

1/2 

stationary  point  at  A  =  r  .  Thus 


g  (x) 


x 


,  ^  1 ,  2  v  2 

1 +  j(x  -r) 


Y 


3/2 


P2<0)  -Y6x2 

e 


h  (x)  =Jl  +  i(x2-r)  2 
h" (A =  r1/2)  =  r 


(4.11) 

(4.12) 

(4.13) 


Therefore 


g  (x)  cox  [h  (x)  t]  dx  = 

4 

0 


3/2 

Y 


p2  (0) 

P^O) 


-  6r 
e 


x 


x  cos ( T 


7T. 

+  t) 


(4.14) 


Thus  for  large  x  eqn.  (4.10)  becomes 

Px(t)  P1(°°)  cos(x+J) 

P1  (0)+P2  (0)"  =P1(0)  +  P2  (0 )  "  C 


where 


' 


, 


i  ii  •  " - 
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(0 ) +p2  (0 ) 


(4.16) 


So  as  time  progresses,  oscillations  of  all  other 
frequencies  will  die  out  faster  than  the  oscillation 
with  frequency  u)=l.  In  other  words  asymptotically 
the  period  of  oscillation  for  any  two  channel  system 
with  any  initial  conditions  is  always  T  =  2tt  . 


For  small  t  the  behavior  of  eqn.  (4.10)  can 


be  quite  different.  As  was  said  in  Section  3.2, 
the  largest  contribution  to  the  time  evolution  equa¬ 
tion  then  comes  from  the  integration  region  where  the 
non-oscillating  factor  of  the  integrand  (i.e.:  g  (x) ) 
has  its  maximum.  Thus  let  us  study  g (x)  and  obtain 
approximations  to  eqn.  (4.10)  for  various  limiting 
cases.  Rewrite  g(x)  as 


g (x)  =  g1  (x)  -  g2  (x) 


(4.17) 


where 


3/2  K2V  ' 


x  2 
—  Y  OX 

e 


(4.18) 


2 


-fix 


2 


x 


e 
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Case  I  p 2  ( 0 )  <<  (0) 

If  p2  (0)  =  0  then  g^  (x)  =  0  and  eqn.  (4.10) 
simplifies  to  eqn.  (1.25).  That  is,  we  are  back  to 
the  one  temperature  system  of  Chapter  One.  Conse¬ 
quently  if  p2(0)  «  p^(0)  we  expect  similar  behavior, 
and  so  we  now  discuss  the  one  temperature  system  and 
the  p2(0)«  Pl(0)  system  together. 

There  are  two  limiting  cases.  If  6  is  large 

then  the  main  contribution  to  the  integral  comes  from 

-1/2 

the  region  around  x  =  6  where  the  factor 

x2e_,Sx2  (4.19) 

attains  its  maximum  value.  Consequently  the  oscilla¬ 
tions  of  eqn.  (4.10)  have  frequency 

W  —  [1  +  -  r>2]1/2  (4-20) 

for  small  times. 

*  2 
—  ox 

On  the  other  hand  if  6  is  small  then  e  ^  1 

for  not  too  large  x  and  the  main  contribution  to  the 

1/2 

integral  is  the  region  around  x  =  r  where  the 
factor 

[1  +  ^-(X2  -  r)2]"1  (4.21) 

attains  its  maximum  value.  Then  the  oscillations  of 
eqn.  (4.10)  have  frequency  w  =  1  for  all  times. 
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Looking  at  Fig.  1.1,  if  we  consider  6=1  large 
and  6  =  0.1  small  then  the  above  arguments  plus 
Kelvin's  stationary  phase  argument  accurately  give 
the  frequencies  for  the  four  examples.  For  instance, 
the  dashed  and  solid  curves  of  the  lower  graph  both 
have  initial  frequencies  uj  =  1  equal  to  the  asymptotic 
frequency.  Thus  these  oscillations  decay  slowly.  The 
dotted  curve  and  the  upper  graph  curve  have  initial 
frequencies  go  *\>  /21  and  go  =  /S'  respectively  which  are 
far  from  the  asymptotic  frequency.  Thus  these  oscil¬ 
lations  decay  quickly.  In  the  upper  graph  we  can  see 
the  asymptotic  frequency  go  =  1  emerging  as  the  frequency 
go  =  /5  decays  away.  Examples  similar  to  these  with 
(0)  /p-^  (0)  /  0  but  small  are  given  in  Fig.  4.1. 

Of  course  if  p2 ( 0 )  >>  p^(0)  then  g2  (x)  is 

negligible  compared  to  g^ (x)  and  we  expect  oscillations 
similar  to  the  above  but  "inverted".  The  above  argu¬ 
ments  for  the  frequencies  are  unchanged  except  that  6 
is  replaced  by  y6.  Examples  of  this  are  given  in 
Fig .  4.2. 


Case  II _ y  large  with  r  =  6  =  1 

To  study  this  case  break  up  eqn.  (4.10)  into 


the  parts 


P-L  (t) 

P1 (0) +p2  (A) 


P1(ot  +  p2('0)  _Il(t)  +  I2(t)  (4-22) 


' 

, 


7b 
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rnu/pr 


8'-  75 


10*700 


Fig.  4 . 1  (a) 


P2  (0) 

p-"(oy  =  0.25;  r  =  5.;  6  =  1 ,  ;  y  =  3 . 

(Compare  this  diagram  to  the  upper  graph  of 
Fig.  1.1.) 


Fig.  4.1(b) 


P2(0) 

p]Jo) 


0.25;  r 


5 .  ;  6  =  1 .  ;  y 


0.5. 


Examples  of  the  two  temperature  model 
with  P2(0)  <<:  P-^(O)- 


— I —  -  —  —  s 


P.HOU  H/C  RHOl  r  0) +-RHO?  C  0  )  )  RHOl  fTl/tRHOl  (  0)+RH0?.(  0  )  )  RHOl  (  T  )/(  RKOl  (  0  J+RH021  01) 

.00  0.Z5  O.SO  0 .75  l.ftO  4). 00  0.2S  0-50  0.7S  l.ftO  cP.*00  0.25  0..S0  0.75 
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Tflu/pr 


(a) 

p1(0) 

PTwr^ioy" 

r  =  5. 

6  =  1. 

Y  =  2. 

(b) 

Px  (0) 

P1  (0) +P2  (0)  " 

r  =  0. 

6  =  0.1. 

Y  =  2. 


(c) 

P1(0) 

Pl(0)+P2(0)=( 

r  =  0. 

6  =  2. 

y  =  0.5. 


Examples  of  the  two  temperature  model  with 
Pl (0)  <<  p2 (0) . 


Fig.  4.2 


■ 


where 


00 


Ii  (t)  = 


(x) cos [h (x) t] dx 


0 


(4.23) 


with  h(x)  given  in  (4.12).  Excluding  1^  (t) ,  for  r  =6  = 
eqn.  (4.22)  is  exactly  the  solid  line  curve  of  the 
lower  graph  of  Fig.  1.1.  By  the  above  arguments  it 
is  a  slowly  damped  oscillation  with  frequency  w2  =  1 
for  all  times.  Thus  we  can  write  1^  (t)  for  limited 
times  as 

I2 (t)  ^  a2  cos  t  (4.24) 

for  some  constant  A2 . 

Looking  at  I, (t) ;  since  y6  is  large,  the  major 
contribution  to  1-^  (t)  comes  from  the  region  around 
x  =  (yS)-1/2.  Thus  I,  (t)  has  frequency 

“l  =  h<i=>  =  [1  +  T<W  '  r)2] 1/2  (4-25) 

/y6  r 

and  we  can  write  I^t)  for  6  =  r  =  1  for  limited  times 
as 

11  o  1/2 

I-  (t)  ^  A-.COS  [1  +  -r  ( — -  “  1)  ]  T  .  (4.26) 

1  l  4  y6 

Since  w  ^  1  this  is  also  a  slowly  damped  oscillation 
but  with  slightly  shorter  period  than  I2 (t) .  Thus 
I  (t)  and  I1(t)  will  interfere  causing  beats.  In  fact 
if  p2  (0)  =  p1(0)  then  ^  A2  and  we  can  use  the  tri¬ 
gonometric  identity 
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cos  a  -  cos 


3=2  sin  2"(a+3) 


sm 


1 


( 3-a) 


(4.27) 


to  write  (4.22)  as 


Px  (t) 

pl (0) +pz  (0) 


% 


p-j^  (°°) 


oo  i  +  oj  2 


roi~p-roy +  A'  sin(-T~)T 


.  ,^1  W2, 

sin  ( — * — )  t 


(4.28) 


For  example  if  y  =  10  then 

sin(1.047T)sin(0.047x) . 

(4.29) 

This  function  describes  Fig.  4.3(a)  very  well.  The 
envelope  period  is 


P-j^  (t) 


Px (0) +p2  (0) 


Px  (°°) 

p^ToTTp^ToT 


+  A 1 


T 

e 


2  7T 

OVOVJ 


^  42 . 5  tt 


(4.30) 


as  the  figure  verifies.  Other  examples  of  this  pheno¬ 
menon  are  shown  in  Figs.  4.3(b)  and  (c) . 
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(a) 


Px  (0) 


P -^  ( 0  )  +  p  2  ( 0 ) 


r  =  1. 
6  =  1. 
y  =  10 


•=0.5 


(b) 


P1(0) 

p^ToT+pJToT 


=  0.6 


r  =  1 
6  =  1 
y  =  5 


TflU/PI 


S.oo 


6.2S 


12.50 


18.75  25.00  31.25  37.50 


TflU/PI 


43.75 


50§00 


(c) 


P1(0) 

P1 (0) +P2  (0) 


•=0.7 


r  =  1. 
6  =  1. 
y  =  10 


Fig.  4.3.  Examples  of  the  two-temperature  model  with 
y  large. 
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CHAPTER  V 


WEAKLY  DEGENERATE  SYSTEM 


§5.1  The  Time  Evolution  Equation 

In  this  chapter  we  study  a  weakly  degenerate 
two-channel  system  of  either  bosons  or  fermions.  We 
start  with  eqn.  (1.23): 


V 


p.  (t)  =  p.  (0) - d3k  i  — tl-  cos(x‘1)-x‘2)  )£]  X 

1  1  <2tt)J  J  2  V2+52 


(5.1) 


where 


✓  2.  2 
.k- 

n(°)(£a))=  [e'^e  1  ±  l]"1 


F-D 

B-E 


.U'-x'2'  -  EJ, 


1/2 


r  _  1  (ed)  _  e<2)) 

2  £k  ek  * 


By  weakly  degenerate  we  mean  that  the  system  approx¬ 
imately  obeys  Maxwell-Boltzmann  statistics.  In 
the  diagram  below  are  plots  of  the  distribution 
functions 


n 


(o) 


(e)  = 


r  a  $e 

[e  e^ 


±  1] 


-1 


[eae3e]  1 


F-D 

B-E 

(5.2) 


M-B 


(where  a=-$y)  as  functions  of  a  (for  fixed  3s). 
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Fig.  5. 


-3  -2  -1  0  +1  +2  +3 


a+8e 


1.  Fermi-Dirac ,  Bose-Einstein ,  and 
Maxwell-Boltzmann  distribution 
functions  as  functions  of  a+$e 
(3e  fixed  arbitrarily  at3e=0). 

Note  that  for  a  large  and  positive 
the  three  distributions  are 
approximately  equal.  Physically 
this  means  that  n^°'  (e)  <<  1  and 
the  effects  of  indistinguishability 
vanish.  (Taken  from  [4]). 


From  Fig.  5.1  we  conclude  that  a  is  large  and 
positive  for  weakly  degenerate  systems.  Using 
the  abbreviation  Z  =  ea  =  e  ^  we  have  that  Z  is 
large.  Thus  we  can  expand  the  distribution  functions 
for  bosons  or  fermions  as  power  series  in  1/Z: 


(o)  ,  \  e 
n  (e)  =  - 


-3e 


-3e  -1 

d±ez  > 

F-D 

B-E 

oo  -n3e 

I  (^)nV 

n=0  Z 

F-D 

B-E 

(5.3) 

Using  (5.3)  we  first  express  p1(0)  as  a  power  series 
in  1/Z. 


P-L  (0)  = 


d3k  n<o)  (e/1  ) 


-e 


2  rru 


(2tt) 


2  e 
k  dk  - 


-n3 


2  2 
fi  k 

2m-. 


0 


I  <TD 

n=0 


n  e 


,n 


F-D 

B-E 


m,  3/2 


=  < - T-> 

2TTBR 


00 

-  y 

7.  L 


(+l)n 


Z  n=0  Zn  (n+1 ) 3/2 


F-D 

B-E 


(5.4) 


Here  we  have  used  the  fact  that 


CO 


0 


,  2  -6x‘ 

dx  x  e 


/? 

4 


1 

T/2 


(5.5) 


Define  the  (dimensionless)  degeneracy  parameter  y: 


.a moJ3\s  -;j: -^e:v  £  >  r  -to!  svxiiaoq 
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y 


P-L  (0) 


,  2  7T  gfi 


2  3/2 
-) 


(5.6) 


Then,  rewriting  (5.4)  in  terms  of  (5.6)  we  get 


(+l)n 


(n+1)  3//2 


F-D 

B-E 


(5.7) 


A  simple  reversion  of  series  [5]  gives  a  power  series 
for  1/Z  in  terms  of  y.  We  get 


+ 


F-D 

B-E 


(5.8) 


Note  that  the  greater  the  degeneracy,  the  larger  the 
value  of  y;  and  also  that  y  increases  with  increasing 
density  and  decreasing  temperature. 

Before  we  can  proceed  we  must  express  the 
distribution  functions  (5.3)  as  power  series  in  y. 
Using  (5.8)  and  simply  Cauchy  multiplying  the  infinite 
series  we  get 


(e) 


1 

Z 


l 

n=0 


(  +  l)n 

„n 


-  (n+1) $£ 
e 


-Be*  2,e  3£  -2$e. 

y  e  ±  y  (  -  e  ) 


3,  -3 Be 
+  y  (e 


/2 


e-23E+(i 


1 

"372  )S  } 


(5.9) 


Using  expression  (5.6)  for  l/p^(0)  and  the  dimension¬ 
less  parameter  substitutions  (1.24)  from  the  two- 
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channel  system  of  Chapter  One  we  get  a  time  evolution 
equation  similar  to  (1.25)  but  with  one  more  parameter - 
the  degeneracy,  y.  That  is 


P1(t) 

p-l  (o  y 


=  1  - 


_2_  63/2 

/? 


oo 


0 


2  . 

x  dx 

,  ^  1 ,  2  77 

1  +  j(x  -r) 


[1  -  cos  t  /l+  x(x2“r)  2] 


r  -fix'  .  /I  -6x‘ 

x  {e  *  y'rv? e 


-2  6x  N 
e  )  + 


(5.10) 


where 


27T$h 


2  3/2 


y  =  pi(0)  (— st') 


T  = 


2V  t 
o 


6  =  —  $V 
m^  o 


(5.11) 


o 


V 


o 


X  = 


2  2 
kZ 

2yV 

o 


U  = 


mlm2 


and 


m2  ml 


1  <  M  Sion.  9fto  /toiw  M Jd  (as  t  )  o^t  7 si  irmra  rco-bfiupe 

* 


§5.2  Example 


An  example  of  a  weakly  degenerate  two-channel 
system  is  shown  in  Fig.  5.2.  It  has  the  same  para¬ 
meters  as  the  solid  line  curve  of  the  lower  graph 
of  Fig.  1.1;  that  is,  r  =  6  =  1.  The  degeneracy 
parameter  is  y  =  0.2. 

The  diagram  shows  that  the  oscillations  for 
the  Fermi-Dirac  particles  are  smaller  and  the  oscil¬ 
lations  for  the  Bose-Einstein  particles  are  larger 
in  amplitude  than  those  for  the  Maxwell-Boltzmann 
particles.  The  physical  explanation  for  this  is  as 
follows . 

Assume  that  the  system  is  at  a  point  in  time 
where  the  density  of  the  a  particles  is  decreasing 
(e.g.  at  t  =  0.5  tt)  and  the  occupation  of  the  states 
of  the  b  particles  is  increasing.  Referring  to 
Fig.  5.1  we  see  that  for  fermions  the  occupation 
number,  say  of  the  b  particles,  is  smaller  for  any 
given  energy  level  due  to  the  Pauli  exclusion  princi¬ 
ple.  Consequently  the  number  of  b  particles  that  can 
be  created  must  be  less,  and  so  the  a  to  b  transfor¬ 
mation  is  inhibited.  For  bosons,  a  larger  number  of 
b  particles  can  be  present  in  states  of  any  energy. 
Consequently  the  transformation  of  particles  a  into  b 


is  enhanced. 
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Fig.  5.2.  r  -  1.;  6  -  1.;  y  -  0.2. ;  Density  oscillations  in  a  weakly  degenerate 
system. 


For  times  when  the  a  particle  density  is 
increasing,  a  similar  argument  shows  that  the  b  to 
a  particle  transformations  again  are  enhanced  for 
bosons  and  inhibited  for  fermions.  Thus  the  density 
oscillations  are  always  larger  for  Bose-Einstein 
particles  and  smaller  for  Fermi-Dirac  particles  than 
for  Maxwell-Boltzmann  particles. 
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CONCLUSION 


We  have  discussed  several  aspects  of  two  and 
three  channel  interacting  systems.  Because  in  our 
model  there  are  no  collisions,  there  is  no 
"Stosszahlansatz "  to  destroy  the  memory  effect  of 
the  system.  Thus  we  get  damped  oscillatory  rather 
than  exponential  approach  to  equilibrium.  We  showed 
that  in  any  n-channel  system  (n  > 1)  we  get  damped 
oscillations  and  that  the  period  of  these  oscillations 
can  increase  greatly  as  the  number  of  channels 
increases.  Thus  this  model  may  explain  some  chemical 
oscillation  reactions  that  have  been  studied  experi¬ 
mentally  . 

We  also  discussed  systems  with  interesting 
initial  conditions,  such  as  the  zero  temperature  limit 
for  fermions  (Chapter  III) ,  two  initial  temperatures 
(Chapter  IV)  and  moderately  high  densities  or  low 
temperatures  (Chapter  V) .  In  each  case  we  discussed 
the  resulting  time  evolution  equation  from  a  mathema¬ 
tical  standpoint,  and  wherever  possible,  gave  physical 
explanations  for  the  behavior  of  the  system. 

There  are  many  aspects  of  this  model  that  should 
still  be  studied.  Some  of  these  are:  Bose-Einstein 
condensation,  the  highly-degenerate  expansion,  and  the 
equilibrium  properties  of  the  model. 
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APPENDIX 


For  reference  purposes,  most  of  the  graphs 
of  the  various  time  evolution  equations  done  to 
this  date  (except  some  of  those  included  in  the 
body  of  the  thesis)  are  assembled  here.  The  inte¬ 
grals  were  solved  numerically  with  an  iterative 

Simpson  rule.  The  n  iteration  partitioned  the 

n 

interval  0  <  x  £  46  2  into  2  subintervals.  If  the 

difference  in  the  values  of  the  integral  computed 

in  the  n^  and  (n-l)^  iteration  was  less  than  10""^ 
t  h 

then  the  n  iteration  value  was  assumed  to  be 
correct . 
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Fio.  A. 1.  r2=  3. ? 

6  =  1. 


m2  m3 


0  .  ;  U1=  1.  ;  U2=  3 .  ; 


0  £  T  £  10  7T 

10  tt  <  t  <  20  71. 


Upper  araph: 
Lowe  r  araph : 
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Fig .  A.  2 .  r0  =  3 .  ;  ru  =  5 .  ;  —  =  —  =  5  •  n  =  i  .  tj  _  -3  . 

2  3  '  m  m  '  ul  '  u2  ~  ' 


6=1. 


Fig.  A. 3.  r2 
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Fig.  A. 4.  r2=  3.; 

6  =  1. 


0 


1. 


TflU/PI 


ml 

Fig.  A. 5.  r2=  3.;  ^  5.;  —  =  ^  =  0.;  U2=  1., 

6  =  0.1. 
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TRU/PI 


Fig.  A. 7.  r2=  3. ; 

6  =  1. 


0. ;  U1=  1. ;  U2=  5. ; 


Note  that  because  U2  is  large,  the  p  oscilla¬ 
tions  damp  out  quicKly. 
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Fig.  A. 10.  r2=  1.; 


Tflu/pr 


6  =  5. 


Tflu/pr 


m,  m. 

Fig.  A. 11.  r0=  1.;  r^=  5.;  —  =  0.;  —  =  0.5. ;  U  =  U  =  6=1. 

2.  3  m2  m3  12 

Note  the  similarity  to  Fig.  2.2(b). 
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Note  the  similarity  to  Fig.  2.6(c).  This  is 
because  U  and  P^,  i  =  1,2,3  are  the  same  at 
x  =  6”"2  for  both. 
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Fig.  A. 14.  r2=  r3= 
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